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ABSTRACT 


The  following  report  documents  a  study  or  low  altitude 
atmospheric  radiation  transmission  conducted  by  a  student-faculty  team 
at  California  State  University,  Northrldge,  School  ot  Engineering. 

The  report  was  generated  under  a  contract  Issued  by  the  Air  Force 
Flight  Test  Center  at  Edwards  Air  Force  Base.  A  simple  computer 
program,  LARTAP,  has  been  developed.  The  simplicity  of  the  program 
maKaS  It  easy  to  modify  to  suit  users  requirements.  Should  detailed 
testing  be  required  to  Improve  program  accuracy,  currently  available 
Instrumentation  has  been  evaluated.  It  has  also  been  suggested  that  a 
tecun i que  Involving  direct  measurement  ot  radiation  transmission  and 
utilization  ot  a  numerical  filtering  technique  may  provide  a  more 
accurate  modei  ot  atmospheric  radiation  transmission  at  a  test  site. 
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INTRODUCTION 


As  radiation  technology  becomes  more  and  more  a  part  of  our 
defense  capabilities,  the  need  arises  to  be  able  to  understand  and 
predict  the  effects  of  the  earth’s  atmosphere  on  radiation 
transmission.  The  need  Is  especially  great  for  the  low  altitude 
regions  of  the  earth's  atmosphere  because  many  systems  are  being 
designed  to  operate  In  this  region.  Particular  examples  of  such 
systems  are  Forward  Looking  Infrared  (FUR)  Lasers,  Low  Light  Level 
Television,  and  Night  Vision  Goggles. 

The  Air  Force  Flight  Test  Center  (AFFT  Edwards  Air  Force 

Base  (Edwards  AFB)  requires  such  knowledge  to  evaluate  sensors  In  the 
visible  through  Infrared  spectrum.  In  order  to  acquire  this  knowledge 
AFFTC  Issued  a  contract  to  the  School  of  Engineering  at  California 
State  University,  Northrldge  (CSUN). 

CSUN  was  to  evaluate  the  transmission  of  radiation  In  the 
atmosphere  of  the  earth.  The  radiation  of  Interest  was  In  the 
wavelenth  region  of  0.35  um  to  14  urn.  The  atmosphere  of  Interest  was 
from  an  altitude  of  0.0  km  (AGL)  to  an  altitude  of  0.7  km  (AGL),  with 
a  maximum  altitude  of  5  km  (ASL).  CSUN  was  also  to  Identify 
atmospheric  components  that  effect  the  transmission  of  radiation  In 
the  aforementioned  wavelength  region. 

The  contract  was  carried  out  at  CSUN  through  the  use  of  a  student 
Design  Clinic.  A  Design  Clinic  Is  an  Independent  study  course  offered 
at  CSUN  where  students  work  under  the  guidance  of  faculty  members  on 
contracTS  Issued  to  the  School  of  Engineering.  The  faculty  and 
student  team  that  worked  on  this  contract  Is  listed  on  the  title  page 


as  authors 


The  project  began  with  a  literature  search.  The  first  goal  of 
the  literature  search  was  to  acquire  and  understand  the  theoretical 
equations  that  describe  the  transmission  of  radiation  In  the  earth's 
atmosphere.  The  second  goal  of  the  literature  search  was  to 
Investigate  the  atmospheric  components  affecting  the  transmission  of 
radiation.  The  final  goal  of  the  literature  search  was  to  evaluate 
presently  existing  computer  software  that  modeled  the  transmission  of 
radiation  In  the  earth's  atmosphere. 

Currently,  the  most  comprehensive  computer  program  that  models 
transmission  of  radiation  In  the  earth's  atmosphere  Is  LOWTRAN. 

LOWTRAN  Is  a  product  of  Air  Force  Geophysics  Laboratory  at  Hanscom 
AFB,  Massachusetts.  This  program  was  first  developed  In  the 
mld-1970's  and  several  modified  versions  have  been  released  since  that 
time. 

Unfortunately,  the  LOWTRAN  computer  code  Is  extremely 
complicated.  This  complication  posed  a  problem  since  AFFTC  had 
expressed  a  need  for  a  computer  program  that  could  be  used  at  a  test 
site.  It  was  determined  that  LOWTRAN  could  not  be  used  to  satisfy 
this  need.  However,  LOWTRAN  contains  an  extremely  comprehensive 
empirical  data  base.  Thus  CSUN  developed  a  new  computer  program, 
LARTAP,  that  Is  much  simpler  than  LOWTRAN,  however  the  new  computer 
program  still  utilizes  the  same  empirical  data  base  as  LOWTRAN. 
Preliminary  test  runs  of  LARTAP  have  shown  agreement  with  LOWTRAN. 
However,  It  Is  suggested  that  further  testing  be  done  to  compare 
LARTAP  to  actual  test  data. 

CSUN  has  also  evaluated  the  AFFTC  test  equipment  to  assess  the 


capability  of  AFFTC  to  measure  atmospheric  components  that  affect 
radiation  transmission.  The  equipment  existing  at  AFFTC  Is  more  than 
adequate  to  support  the  current  version  of  LARTAP.  However,  should 
the  need  arise  for  more  comprehensive  measurements,  suggestions  are 
made  for  possible  acquisitions. 

Finally,  suggestions  are  made  for  possible  future  work.  As  was 
mentioned  previously,  more  testing  needs  to  be  performed  to  verify  the 
LARTAP  computer  program.  One  of  the  short  comings  of  both  LARTAP  and 
LOWTRAN  Is  the  Inability  to  handle  the  possible  affects  of  solar 
radiation.  It  Is  suggested  that  a  model  be  developed  to  account  for 
this  affect.  The  resulting  model  could  easily  be  Incorporated  In 
LARTAP. 


SECTION  I 


I 

BASIC  CONCEPTS  IN  RADIATION  TRANSMISSION 

If  a  source  emits  radiation  Into  a  medium  (such  as  the  atmosphere 
of  the  earth)  then  the  Intensity  of  the  emitted  radiation  will  be 
changed  by  the  medium.  The  situation  Is  presented  qualitatively  In 
figure  t.  The  radiative  Intensity  from  the  source  of  Interest  Is 
changed  through  four  mechanisms.  Radiation  may  be  absorbed,  radiation 
may  be  scattered  out  of  the  line  of  sight,  radiation  may  be  emitted 
and  radiation  from  sources  other  than  the  one  of  Interest  may  be 
scattered  Into  the  line  of  sight. 

It  Is  demonstrated  In  reference  1  that  the  change  In  radiative 
intensity  of  wavelength  along  some  path  S  through  a  medium  Is  given 
by  the  following  differential  equation. 

d  I  -  d  I;,  -  a  X  -  cnI  + 

d  5  A-rr 

where : 

S  -  path 

I  -  monochromatic  radiative  Intensity  at  S 

1^  -  blackbody  Intensity  using  the  temperature  of  the  medium  at  S 
.  and  wavelength 

a  -  absorption  coefficient  of  the  medium  at  S  and 

<t\  -  scattering  coefficient  of  the  medium  at  S  and 

w  -  sol  Id  angle 

1  -  phase  function  of  scattering 

I*  -  radiative  Intensity  along  S  from  sources  other  than  the  one 
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of  Interest. 

The  first  term  to  the  right  of  the  equal  sign  accounts  for 
emission  by  the  medium.  Note  that  Klrkoff's  Law  has  been  Invoked  In 
the  derivation  of  equation  (1)  since  the  absorption  coefficient  Is 
used  In  the  emission  term.  The  second  and  third  terms  account  for 
attenuation  due  to  the  absorption  and  scattering  out  of  the  line  of  • 
sight,  respectively.  The  fourth  term  accounts  for  the  scattering  of 
radiation  from  sources  other  than  the  one  of  Interest  Into  the  line  of 
sight. 

The  scattering  coefficient  and  phase  function  In  the  fourth  term 
may  be  Interpreted  qualitatively  as  follows.  The  scattering 
coefficient  determines  the  total  amount  of  radiation  from  sources 
other  than  the  one  of  Interest  that  Is  scattered.  The  phase  function 
accounts  for  the  proportion  of  radiation  that  Is  scattered  Into  the 
line  of  sight.  A  good  example  of  I*  source  of  radiation  would  be 
solar  Intensity. 

Unfortunately,  In  order  to  properly  evaluate  the  last  term  In 
equation  (1)  It  would  be  necessary  to  know  the  local  intensity  along 
the  line  of  sight  from  all  possible  radiative  sources  around  and  In 
the  region  of  the  medium  under  consideration.  This  Is  quite  a 
formidable  task.  Briefly  It  Involves  transforming  equation  (1)  Into 
an  I ntegro-part I  a  I -d I f ferent I al  equation.  Then  this  equation  Is 
solved  for  an  Intensity  field  function  that  may  be  evaluated  at  any 
point  In  the  region  of  the  medium  under  consideration.  In  general 
this  I ntegro-part I al -d I f ferent I  a  I  equation  may  be  solved  only  by  a 
numerical  technique  such  as  finite  element. 

Thus  In  the  Interest  of  simplicity  the  last  term  of  equation  (1) 
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was  dropped.  It  should  be  pointed  out,  however,  that  for  a  I  I ne  of 
sight  that  extends  from  the  surface  of  the  earth  upwards  I*  source 
term  resulting  from  the  scattering  of  solar  radiation  may  be  quite 
significant.  It  Is  felt  that  a  future  effort  should  be  undertaken  to 
find  a  simplified  method  to  model  this  effect. 

Dropping  the  last  term,  equation  (1)  may  be  written  as: 


d_X  +  6  I  =  alb 
d  S 

where: 

€  *  +  a  -  extinction  coefficient 

Equation  (2)  Is  a  linear,  ordinary  differential  equation  which 


has  the  following  solution. 


IL  =  I. 


-(  6(u)du  •  f  "Ce  ^  u 

^  \  IvCOoCs*)  ds* 


(.3 


where  I,,  and  I0  are  the  following  Dlrlchlet  boundary  conditions. 

I0  -  monochromatic  Intertsl+y  at  the  source  of  Interest, 

IL  -  monochromatic  Intensity  at  a  point  a  distance  L  along  S. 

Equation  (3)  may  be  Interpreted  qualitative  as  follows.  Consider 
the  case  where  the  medium  Is  cold  such  that  0,  then  equation  (3) 

would  reduce  to: 


iL  =  i. 


-  C  €.  (  u)  d  U 

Jo 


Thus  the  exponential  with  an  Integral  augument  represents  the 


attenuation  of  the  Intensity  I0  along  the  line  of  sight  extending  from 
0  to  L.  This  exponential  with  an  Integral  augument  defines  the 
transmittance  of  the  path  from  0  to  L.  Note  that  the  transmittance 


will  vary  from  zero  to  one. 

The  second  term  Is  a  little  more  compt Icated.  Consider  a 
differential  path  length  at  S*  some  where  along  S  between  0  and  L. 
Over  this  differential  path  length  dS*  radiation  will  be  emitted  by 
the  medium.  The  exact  amount  Is  given  by  (S*)a(S*)dS*.  However, 
this  amount  will  be  attenuated  between  S*  and  L  where  the 


transmittance  Is  given  by: 


r  L 

-  ^  e  (  u)  d  u 


Summing  all  the  contributions  from  every  dS*  between  0  and  L  gives  the 
second  term  In  equation  (3). 

Equation  (3)  may  be  written  as  follows: 


-C  £  Cu)  cl  u 

Iv.  =  Io  «T° 

+  Wes*)  e 

J  -  ^  o(w)  d  ^ 
e  i 

Define  Ta  and  T«  as  follows: 

-  ^  €  (  u)  d  u 

t«  (x.  y)  =  e 


-  V  a(u)  du 


-  ^  a  (u)  d  ^ 


T«  (*,  'i')  -  G 


Then  equation  (4)  may  be  written  as: 

II  «IoT.(o,l)  +■  f  Iw(s«)  TeC&lL)  d  T(s»,  L)  ds* 

Jo  To  (s*.  l)  d  s* 


At  first  glance  It  may  appear  that  equation  (5)  Is  more 


"  .'*.**.  I  •  *  •  ’  •  *  • 
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complicated  than  equation  (3),  however,  the  numerical  analogue  of 
equation  (5)  will  prove  to  be  very  convenient  for  evaluating  Iu  with 
a  computer  program.  The  numerical  analogue  of  equation  (5)  is: 

Iu  =  l„  rt(o,  O 


ItCv.,)  +  i>CO  ZiULO 

l)  t0(s“,l) 


x  T.  (S**,  ,  L)  -  T.Cs'.L) 


where  the  path  from  0  to  L  has  been  divided  Into  n-1  distances  of 

(  $  i+i  “  -*  i.  )  • 


SECTION  II 

RADIATION  PROPERTIES  OF  GASES  AND  PARTICLES 

In  section  I  the  concepts  of  absorption,  emission  and  scattering 
were  Introduced.  The  purpose  of  this  section  Is  to  discuss  these 
concepts  In  more  detail.  However,  before  discussing  the  microscopic 
effects  producing  these  three  phenomenon  It  Is  helpful  to  review  the 
continuous  matter  assumption. 

Equation  (1)  assumes  continuity  of  matter  In  the  following  sense. 
A  macroscopic  differential  volume  of  matter  contains  an  Infinite 

number  of  molecules;  such  that.  In  equilibrium  although  the 

% 

microscopic  states  of  the  molecules  are  continuously  changing,  the 
macroscopic  state  of  the  matter  In  dV  Is  the  time  average  of  the 
continuously  changing  microscopic  states  of  the  molecules  within  dV. 

In  order  to  understand  absorption,  emission  and  scattering  It  Is 
necessary  to  utilize  Quantum  Mechanics.  From  a  Quantum  Mechanical 
point  of  view  the  radiative  transfer  problem  becomes  one  of 
determining  the  characteristics  of  photons.  In  particular  the  local 
number  density,  energy  distribution  and  direction  of  motion.  If  these 
three  characteristics  are  known  then  the  macroscopic  quantity. 
Intensity,  may  be  determined.  Any  or  all  of  these  characteristics  are 
altered  by  Interactions  between  the  photons  and  the  medium  components 
(atoms,  molecules  and  particulates). 

Absorption  occurs  when  a  photon  encounters  a  component  of  the 
medium  and  Is  absorbed  thereby  Increasing  the  energy  of  that 
component.  Emission  occurs  when  a  component  emits  a  photon  thereby 
decreasing  the  energy  of  that  component.  Scattering  occurs  when  a 
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photon  encounters  a  component  of  the  medium  and  the  direction  of 
motion  of  the  photon  Is  changed. 

According  to  Quantum  Mechanics  a  single  medtum  component  may 
exist  In  any  one  of  a  number  of  accessible  energy  states.  This 
concept  Is  displayed  In  figure  2.  The  ground  state  Is  the  lowest 
energy  level  that  a  medium  component  may  have.  The  ionization  state 
Is  the  energy  required  to  Ionize  the  medtum  component.  In  figure  2, 
two  other  accessible  energy  states  are  shown  between  ground  and 
Ionization.  For  an  actual  medium  component  there  may  be  any  number  of 
accessible  energy  states  between  ground  and  Ionization. 

Consider  a  single  gas  molecule  that  has  accessible  energy  states 
as  shown  In  figure  2.  According  to  Quantum  Mechanics  the  energy  of  a 
photon  Is  hv*  where  h  !.s  Planck's  Constant  and  T  Is  frequency.  Thus  In 
order  for  this  gas  molecule  to  move  from  energy  state  E,  to  energy 
state  It  must  absorb  a  photon  of  the  following  frequency. 

V  =  E.  -  E, 

h 

If  there  were  no  Interaction  between'  mo  I  ecu  I es  and  no  other 
effects  the  absorption  c6eff!c!ent  for  a  large  number  of  these  gas 
molecules  would  be  a  single  vertical  line  at  frequency  V.  The 
magnitude  of  the  absorption  coefficient  would  depend  of  the  number  of 
gas  molecules  within  the  medium  at  energy  .level  E( .  However,  In  the 
case  of  a  real  gases  collisions  between  molecules  cause  the  accessible 
states  displayed  In  figure  2  to  shift.  Thus  for  a  real  gas  the 
absorption  coefficient  Is  a  distribution  around  frequency  v\  The 
spread  of  the  distribution  depends  on  the  number  of  molecules 
experiencing  a  collision.  The  number  of  gas  molecules  experiencing  a 
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collision  Is  related  to  pressure.  As  was  previously  mentioned  the 
magnitude  of  the  absorption  coefficient  at  frequency  v  depends  on  the 
number  of  molecules  In  energy  state  E,.  This  number  Is  dependent  on 
the  density  and  temperature  of  the  gas. 

Thus  the  following  conclusions  can  be  drawn  about  the  absorption 
coefficient  of  a  gas.  First,  the  absorption  coefficient  will  depend 
on  the  local  thermodynamic  state  of  the  gas.  Second,  the  absorption 
coefficient  for  a  gas  Is  a  strong  function  of  Incident  radiation 
frequency.  The  principles  presented  thus  far  are  displayed  In  figure 
3. 

Finally,  If  the  schematic  displayed  In  figure  2  were  for  a  real 
gas,  then  a  plot  of  absorption  coefficient  versus  all  frequencies 
would  appear  as  three  distributions  of  th*  type  displayed  In  figure  3. 
These  distributions  would  be  centered  at  the  following  frequencies: 


b  b 


The  relative  magnitude  of  the  absorption  coefficient  at  each  of  these 
frequencies  would  depend  of  the  thermodynamic  state  of  the  gas. 

Attention  Is  now  turned  to  the  medium  components  that  are 
referred  to  as  particulates  or  aerosols.  These  components  are  solid 
particles  that  are  suspended  In  the  gaseous  component  of  the  medium. 
Examples  of  particulates  are  dust  particles,  water  droplets,  etc. 

It  Is  Important  to  note  at  this  point  the  size  of  particulates 
and  how  It  relates  to  the  continuum  assumption  discussed  at  the 


90 


beginning  of  this  section.  In  order  to  treat  the  particulates  using  a 
continuum  equation  such  as  equation  (1),  It  must  be  assumed  that  the 
dimensions  of  the  problem  under  consideration  are  large  enough  such 
that  a  differential  volume  dV  contatns  a  large  number  of  particulates. 

Where  as  a  gas  exhibits  an  absorption  coefficient  that  Is  a 
strong  function  of  frequency,  particulates  tend  to  exhibit  absorption 
coefficients  that  are  somewhat  smoother  functions  of  frequency.  The 
reason  Is  that  a  particulate  Is  composed  of  a  number  of  closely 
Interacting  molecules.  Figure  2  for  such  particulates  becomes  a 
continuum  of  accessible  states.  Thus  the  molecular  collision-pressure 
relationship  that  Is  significant  for  gas  has  essentially  been  carried 
to  an  extreme  In  the  case  of  particulates  resulting  In  a  smoothing  of 
the  frequency  dependence. 

The  absorption  coefficient  for  particulates  also  depends  on  the 
thermodynamic  state  of  each  particle  since  this  Is  related  to  the 
energy  state  of  each  molecule  comprising  the  particle.  However  this 
dependence  Is  negligible. 

Thus  the  absorption  coefficient  for  particulates  In  the  medium  Is 
a  function  of  frequency  and  particle  density.  The  density  dependence 
Is  expanded  to  Include  particulate  size  distribution,  since  larger 
particles  will  tend  to  absorb  more  radiation  than  a  smaller  particle. 
This  Is  due  simply  to  the  larger  cross-section  presented  to  the 
Incident  radiation  by  larger  particles. 

If  a  photon  encounters  a  molecule  and  the  energy  of  the  photon  is 
such  that  the  molecule  cannot  absorb  It,  then  that  photon  will  be 
scattered  Into  a  different  direction  as  a  result  of  the  collision. 

In  the  case  of  particulates  the  scattering  phenomenon  Is 
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complicated  by  the  fact  that  when  a  photon  encounters  a  particulate  It 
encounters  the  surface  molecules  of  the  structure.  Thus  It  Is  these 
molecules  that  determine  the  resulting  motion  of  the  photon.  For 
example,  a  photon  may  encounter  a  surface  molecule  that  scatters  the 
photon  toward  a  molecule  on  the  same  particulate  that  finally  absorbs 
this  photon. 

As  can  be  seen  by  this  last  example  the  mechanism  that  govern 
scattering  can  get  quite  complicated.  However,  scattering 
coefficients  of  molecules  and  part  I cu I ates  exh I b 1 1  the  same  dependence 
characteristics  as  the  absorption  coefficients. 


SECTION  III 


RADIATION  MODEL 


Given  the  tools  developed  In  sections  I  &  II  It  Is  now  possible 
to  formulate  a  mathematical  model  to  analyze  the  transmission  of 
thermal  radiation  through  the  earth’s  atmosphere.  Such  a  model  has 
been  developed  by  C.S.U.N.  In  the  form  of  a  computer  program.  This 
section  documents  this  computer  program. 

The  model  was  designed  tc  function  under  the  following 
constraints.  The  radiation  to  be  considered  by  the  model  Is  In  the 
wavelength  range  from  0.25  um  to  25  urn.  The  region  of  the  atmosphere 
to  be  considered  extends  from  ground  level  to  an  altitude  of  0.7  km 
above  ground  level  with  a  maximum  of  5  km  above  sea  level. 

The  objective  of  the  computer  program  Is  to  evaluate  equation  (6) 
of  section  I.  From  the  Information  presented  In  section  II  It  can  be 
concluded  that  the  absorption  and  extinction  coefficients  will  be  a 
function  of  the  local  thermodynamic  state  of  the  atmosphere  and  the 
frequency  of  the  Incident  radiation.  In  general  the  thermodynamic 
state  of  the  atmosphere  will  vary  In  all  directions,  however,  for  this 
report  It  Is  assumed  only  to  be  a  function  of  altitude  measured 
vertically  to  the  earth. 

The  atmosphere  Is  composed  of  a  large  number  of  components.  It 
Is  assumed  that  the  total  absorption  coefficient  of  the  atmosphere  Is 
the  summation  of  the  absorption  coefficients  of  each  of  the 
components:  n 
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This  assumption  Is  also  made  for  the  total  scattering 
coef f iclent: 


o\ 


tot 


This  Implies  that  the  total  extl  nctlon'  coef  f  Iclent  as  defined  In 
section  I  Is: 


tot 


These  assumptions  lead  to  the  following  conclusions  concerning  T 
and  Te  as  defined  In  section  I. 


Also,  based  on  the  above  assumptions  and  definitions  In  section 
I,  To-  and  for  atmospheric  component  I  are  functions  of  path 
length,  the  thermodynamic  state  of  component  I  along  that  path  length 
and  the  frequency  of  the  Incident  radiation. 

Thus  If  empirical  data  were  obtained  to  establish  functional 
relationships  for  To^  and  T*-,  then  these  relationships  and  equation 
(6)  would  provide  the  desired  result.  This  Is  exactly  the  approach 
that  was  taken  In  formulating  the  model  developed  at  C.S.U.N. 

Fortunately,  the  required  empirical  data  has  already  been 
obtained  for  use  In  a  computer  program  known  LOWTRAN.  This  computer 
program  was  developed  by  the  Geophysics  Laboratory  In  Cambridge,  Mass. 
The  objective  of  LOWTRAN  Is  to  solve  problems  similar  to  the  one 


considered  in  this  report.  The  obvious  question  now  arises  as  to  why 
It  was  necessary  to  develop  a  new  program? 


The  answer  to  this  question  lies  In  the  fact  that  LOWTRAN  was 
developed  to  handle  problems  more  general  than  the  one  under 
consideration  In  this  report.  The  consequence  of  this  generality  Is 
complication.  Thus  by  developing  a  new  model  and  extracting  the 
pertinent  data  from  LOWTRAN  as  required,  C.S.U.N.  was  able  to  develop 
a  simpler  computer  program  that  was  geared  only  to  the  problem  at 
hand. 

According  to  LOWTRAN  documentation  (reference  2),  rYa(x,y)  and 
Te  L  (x,y)  may  be  expressed  as  follows: 


£u v\  c  o 


Cr;  (r) 


UJL  (>)  <A  $ 


where : 

T  -  frequency 

Crc~  spectral  coefficient  for  component  I 
h  -  altitude  measured  vertically  from  the  ground 
s  -  path  (x  to  y) 

U  i  ~  absorber  or  scattering  amount  of  component  I 

There  Is  an  assumption  that  needs  to  be  made  In  order  to  utilize 
the  empirical  data  from  LOWTRAN.  The  assumption  is  related  to  the 
fact  that  equation  (6)  from  section  I  applies  only  for  radiation  of  a 
discrete  frequency.  Unfortunately,  It  is  virtually  Impossible  to 
obtain  the  required  empirical  data  at  discrete  frequencies.  The  data 
obtained  for  LOWTRAN  was  acquired  by  Instruments  that  essentially 


obtained  the  average  values  of  /Va  and  over  a  narrow  band  of 
frequencies.  Thus  in  order  to  utilize  equation  (6)  the  assumption 
must  be  made  that  this  equation  Is  also  valid  over  a  narrow  band  of 
f requenc les . 

LOWTRAN  provided  data  for  the  following  atmospheric  components: 

Ha0,  0*,  C0Z,  CH+,  N*0,  CO,  0-w  Na 
and  also  five  general  categories  of  particulates  (aerosols) 

1.  Maritime 

2.  Urban 

3.  Rural 

4.  Advectlon  Fog 

5.  Radiation  Fog 

The  components  CO* ,  CH*,  Nx0,  CO,  and  0*.  were  lumped  as  one 
component  which  Is  referred  to  as  the  uniformly  mixed  gases.  The 
absorber  amount  for  this  component  Is  given  by: 

O.  IS 

U  =  p  a-?v  is  P 
1013.0  T  _\CUS.O 

where : 

P  -  local  total  pressure  (mb) 

T  -  local  temperature  (K°) 

For  Incident  radiation  having  a  wavelength  In  the  region  of  4  urn, 
Nx  molecules  experiencing  a  collision  provide  a  continuum  of 
accessible  energy  states  which  are  capable  of  absorbing  radiation  In 
this  region.  This  effect  Is  considered  by  the  program  and  the 
absorber  amount  Is  given  by: 


where: 

P  -  local  total  pressure  (mb) 
T  -  local  pressure  (K°) 


The  atmospheric  component  Ha0  Is  treated  as  three  separate 
components  In  order  to  artificially  account  for  three  separate 
absorption  characteristics.  The  first  characteristic  Involves  the 
standard  type  of  gaseous  absorption  that  was  discussed  In  section  ||. 
The  absorber  amount  for  this  component  Is  given  by: 


where: 

P  -  local  total  pressure  (mb) 

T  -  local  temperature  (K°) 

UJHZ  -  local  density  of  Ha0  (gm/ms) 


The  other  two  absorption  characteristics  attributed  to  H^O 
Involve  continuum  absorption.  The  two  continuums  are  In  the 
wavelength  regions  of  3.5  urn  to  4.2  urn  and  8  urn  to  11  urn. 
Unfortunately,  theoretical  explanations  for  the  existence  of  these  two 
regions  of  continuum  absorption  does  not  presently  exist.  However 
they  are  accounted  for,  and  the  absorber  amounts  are  calculated  as 
fol lows: 
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8.0  wn  to  11.0  um 


where : 

PH  0  -  local  partial  pressure  of  H40  (atm.) 
P  -  local  total  pressure  (atm.) 


T  -  local  temperature  (  K) 

IjJW2  -  local  density  of  H40  gm/cm  -  km 
3.5  urn  to  4.2  urn 


u  =  u 


HI  PH, 


xo  +  0.I2.(P  -  P*l0)  EXP  A<?A-  - 


where  nomenclature  and  units  are  the  same  as  8.0  um  to  11.0  um  region. 


The  atmospheric  component  09  Is  treated  as  two  separate 
components  In  order  to  artificially  account  for  two  separate 
absorption  characteristics.  The  first  0,  component  accounts  for  the 
standard  type  of  gaseous  absorption  that  was  discussed  In  section  II. 
This  component  will  be  referred  to  as  Ozone  I.R.  Where  the  I.R. 
refers  to  Infrared.  The  absorber  amount  Is  given  as  follows: 


where : 


P  -  local  total  pressure  (mb) 
T  -  local  temperature  (  K) 

Uot  -  density  of  0?  (gm/m  ) 
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The  second  03  component  accounts  for  the  absorption  of  photons  by 
the  electronic  structure  of  the  03  molecule.  This  type  of  absorption 
was  not  discussed  In  section  II  because  It  would  have  required  a 
complete  dissertation  In  Quantum  Mechanics.  Briefly,  absorption  at 
the  electronic  level  of  molecules  Is  associated  with  high  frequency 
Incident  radiation  which  Is  referred  to  as  ultraviolet  (U.V.)  (1  um  to 
0.01  um).  The  absorber  amount  Is  given  as  follows: 

U  *  46.6667  UJol 
where : 

bJ0l-  density  of  03  (gm/m^) 


According,  to  the  LOWTRAN  documentation,  the  total  effect  of 
scattering  by  the  molecular  atmospheric  components  discussed  thus  far 
can  be  lumped  together.  The  scattering  amount  Is  given  as  follows: 


where : 
P 
T 


u  =  P  a. 73.  is 
10  13.0  T 

local  total  pressure  (mb) 
local  temperature  (K°) 


The  spectral  coefficients  for  each  of  the  attenuation  components 
discussed  thus  far  are  presented  In  figure  4  to  figure  11.  In  these 
figures  the  Independent  spectral  quantity  Is  the  wave  number  (*])  and 
It  Is  defined  as  follows: 

kj  =.  i .  o/  A 

where  Is  the  wavelength.  Thus  a  wave  number  of  rj  cm  '  Is  equivalent 
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FIGURE  4 

UNIFORMLY  MIXED  GASES 


FIGURE  5 

NITROGEN  CONTINUUM  (4  UM) 


FIGURE  6 
HSO  GAS 
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FIGURE  7 

H20  CONTINUUM  (8  UN  TO  11  UN) 


to  a  wavelength  of: 


r>)  -  I  O  *  LA  m 

•7 

It  Is  worthwhile  to  note  that  all  -of  the  thermodynamic  properties 
required  to  evaluate  each  WL  are  functions  of  h.  Thus  each  W;_  Is  a 
function  of  h. 

As  was  mentioned  earlier,  LOWTRAN  provided  empirical  data  for 
five  general  categories  of  particulates  (aerosols).  They  are  as 
f ol I ows : 

1 .  Maritime 

2.  Urban 

3.  Rural 

4.  Advect I  on  Fog 

5.  Radiation  Fog 

The  Marttfne  category  applies  to  oceanic  regions  and  coast  line 
areas.  The  particulates  result  from  sea  salt  particles  which  are 
produced  by  the  evaporation  of  sea  spray  droplets. 

The  Urban  category  applies  to  Industrial  and  urban  areas  where 
the  majority  of  particulates  result  from  combustion  processes. 

The  Rural  category  Is  Intended  to  represent  the  aerosol 
conditions  found  In  continental  areas  which  are  not  directly 
Influenced  by  urban  partlcu!*  '<  sources. 

When  the  air  becomes  nearh'  laturated  with  water  vapor,  fog  can 
form.  Saturation  of  the  air  can  occur  as  a  result  of  two  different 
processes;  the  mixing  of  air  masses  with  different  temperatures 
(advectlon  fog)  or  by  cooling  of  the  air  to  the  point  where  the 
temperature  approaches  the  dew  point  (radiation  fog).  It  should  be 
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mentioned  that  a  cloud  cover  Is,  nothing  more  than  a  fog  that  does  not 
contact  the  ground. 

Based  on  the  Information  presented  In  section  II,  the  extinction 
and  absorption  coefficients  for  particulates  are  a  function  of 
particulate  concentrat I  on,  s I ze  distribution  and  frequency  of  incident 
radiation.  The  spectral  and  size  distribution  data  was  obtained  from 
LOWTRAN  for  the  five  categories  of  particulates.  This  data  is 
presented  In  figure  15  through  figure  24.  The  particle  size 
distribution  for  each  of  the  five  categories  of  particulates  was 
obtained  empirically.  It  only  remains  to  account  for  variations  In 
particulate  concentration.  Note  that  In  each  case  the  data  has  been 
normalized  such  that  the  extinction  coefficient  Is  one  at  a  wavelength 
of  0.55  um.  The  reason  Is  to  Incorporate  the  particulate 
concentration  utilizing  a  quantity  called  the  Metero I og I ca I  Range  (V). 
The  quantity  V  Is  the  distance  required  for  the  radiant  Intensity  of 
some  source  to  drop  to  2 %  of  the  original  source  Intensity  at  a 
wavelength  of  0.55  wm.  If  the  source  Intensity  Is  such  that  emission 
by  the  particulates  Is  negllgable,  then 

-  ve 

%  (o,  \l)  -  0.02.  =  <2. 

or 

&  =  3  ,  q  i  2 

V 

where 

B  -  particulate  extinction  coefficient  at  0.55  urn  (km  1 ). 

The  obvious  question  arises  as  to  the  effect  of  other  atmospheric 


38 


components  on  B.  The  answer  lies  In  the  specified  wavelength  of 
0.55  um,  which  Is  approximately  equal  to  a  wave  number  of  18000  cm"1 . 

By  reviewing  figure  4  through  figure  11  It  can  be  seen  that  the 
only  atmospheric  components  affecting  T6(0,V)  at  0.55  otm  are  03  and 
molecular  scattering.  The  concentration  of  03  at  low  altitudes  Is 
virtually  negligible.  It  was  only  considered  In  this  model  for 
completeness.  Also  the  effect  of  molecular  scattering  at  0.55  urn  Is 
negligible.  Thus  the  quantity  V  depends  on  the  particulate 
concentrat I  on . 

The  reason  for  normalizing  the  particulate  data  to  an  extinction 
coefficient  of  one  at  0.55  um  should  now  be  obvious.  By  measuring  or 
estimating  a  value  of  V  and  selecting  the  appropriate  particulate 
category,  the  extinction  and  absorption  coefficients  at  any  other 
wavelength  may  be  calculated  as  follows: 

flL  =  3.  <3  I  2  C  \r  a  -L  ^Hz) 

V 

£  i  =•  3.412,  £  «*e  C  ,  UJ^z.) 

V 

where : 

(3  -  wavelength 

Whz  -  h^O  density 

V  -  metero I og I ca I  range 

C^Q  and  C„6  are  presented  In  figure  15  through  figure  24. 

Note  that  C y*  and  C[r6  are  functions  of  the  H20  density.  This  Is 
shown  In  the  figures  by  relative  humidity.  The  reason  for  this 
functional  relationship  Is  that  the  H20  molecules  tend  to  accrete  to 
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FIGURE  23 
RADIATION  FOG 
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RADIATION  FOG 
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the  surface  of  the  particulates.  Thus  the  Incident  photons  encounter 
H^O  molecules.  The  higher  the  more  accretion  takes  place. 

It  Is  worth  mentioning  some  more  Information  about  the  fog 
models.  For  thick  fogs  (clouds)  where  V  Is  less  than  200m  the 
extinction  coefficient  Is  virtually  Independent  of  wavelength.  These 
conditions  are  modeled  best  by  the  advectlon  fog.  For  moderate  fogs 
200  m  <  V  <  1000  m  the  radiation  fog  should  be  used.  For  thin  fogs, 

1  km  <  V  <  2  km  a  marltlne,  urban  or  rural  particulate  model  should  be 
used  with  a  99$  relative  humidity. 

Utilizing  equation  4b  from  section  I  and  the  above  equations  for 
the  particulate  absorption  and  extinction  coefficient,  To,(x,y)  and 
Tfe(x,y) f°r  a  particulate  atmospheric  component  may  be  expressed  as 
fol lows: 


where 

(>)  -  3.  q  i  a. 

v 

'a] hi  “  shows  the  dependence  of  the  spectral  coefficient  on  the  HA0 
density  (relative  humidity) 
it  -  frequency 

(jJhz)  “  from  figure  14  through  figure  24.  (extinction  and 
absorpt 1  on ) 

In  -  altitude 

The  spectral  coefficient  for  a  particulate  atmospheric  component 
Is  part  of  the  Integrand  due  to  Its  dependence  on  Hz0  density  which  Is 
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dependent  of  h.  This  Is  In  contrast  to  the  equivalent  expression  for 
molecular  atmospheric  components  (see  equation  8).  The  above  equation 
can  be  placed  In  the  same  form  as  equation  8  by  evaluating  the 
spectral  coefficient  at  the  average  Ha0  density  between  x  and  y.  Then 
Ta ■(x.ij')  and  T€-(x,y)  for  the  particulate  atmospheric  component  may  be 
expressed  as  follows: 

_  _  r  ^ 

•Pw^chion  (u-  ,  UJ  hz.  )  \  WtWJs  . 

Jx  J 


where 


^  HZ.  —  Z.  (  x)  ■+■  ^  H  Z.  C  4  ) 


Based  on  the  data  presented  thus  far  In  this  section  It  Is  now 
possible  to  evaluate  the  following  term  for  each  of  the  atmospheric 


components 


!<c  =  C  vr c  C  O 


U  L  C  ^0  d  s 


All  that  remains  Is  to  know  the  functional  relationship  between 
T6  •  (x,y),  'Tai(x,y),  and  K;.. 


For  the  Unlformlly  Mixed  Gases 

T€=T*  =  EXP  [  - 

For  the  Nitrogen  Continuum 


0.55  5  " 

O.  04H  K 

_ 


For  the  H^O 


Te  =  Tfl  =  EXP 


T*  =  T«  =  EXP 


-  K 


o . 06  K 


O.  555 


For  the  H,0  Continuum  (8.0  urn  to  11.0  urn) 


For  the  H^O  Continuum  (3.5  um  to  4.2  um) 

Te  =  T.  =  E  X  P  [-  «  ]  (, 

For  the  0^  I .R.  Region 

O.  &  i  (a 

Tfe  =  To  =  E  X  P  [-  O.  0  64-3  K  J  Os) 

For  the  0a  U.V.  RegJon 

r6  =  t«  =  exp  [-  k]  CO 

For  Molecular  Scattering 

=  EXP  [-  K  ]  0  7) 

Ta  =  0.  O  (l  8) 

For  Particulates 

%  EXP  [-K6]  ■  0*0 

T«  =  EXp['Ka]  (ao) 


Note  that  for  a  particulate  component  there  are  two  values  for  K 
where  KQIs  obtained  by  evaluating  equation  (9)  for  absorption  and  K6 
Is  obtained  by  evaluating  equation  (9)  for  extinction. 

Consider  the  situation  depicted  In  figure  25,  where 
h,  -  altitude  of  point  x 
hA  -  altitude  of  point  y 
T(h,)  -  temperature  at  h, 

T(ha)  -  temperature  at  ha 
UJ C ( h i  )  -  absorber  amount  of  component  l  at  h, 

CdtOij.)  “  absorber  amount  of  component  I  at  ha 

thus 

c\  5  x  d  In 

COS  Q 


CO  \  Wi.(h)ds  =  CO  \  l^c(h)  d  In 

Cos  &  j v 


assume  that  U{.(h)  Is  a  semi-log  function  of  altitude,  then 

h  -  \n . 


U,Ch)  =  UJ-CK) 


u  i  ( O 


<1*-  Vi, 


c, 


Lwt  Cw.)J 

(0  ^  ^  CO  d  ^  =  [uc ( hx)  -  t»j ^  Ch,)]  [  y  -  *]  c a--  (it) 


°  3* 


as  h,  -*  hx 


UJl  C^x) 

Wc  (k)  J 


Co-^ir)  (  UJC  C^)  cl'0  - ^  l*J(.  Ck)  "  XJ  Co\  (y) 


from  the  last  two  expressions  and  equations  (10)  through  (20)  It  Is 
possible  to  evaluate  equations  (7a)  and  (7b)  for  TeT0T(*»y)  ancl 
T„tot  (x,y).  Thus  equation  (6)  of  section  I  can  now  be  evaluated  where 
T(h,  )  and  T(ha  )  would  be  used  In  the  black  body  function.  The 
quantity  n  for  the  numerical  analogue  of  the  emission  Integral  In 
equation  (6)  would  be  2  for  this  particular  example,  and  also  note 
that  In  general . 

T€toT  C  4  >  y)  =  '*° 

tot  C y  ,  y  )  l  -  O 

The  above  has  been  a  simplified  example  of  how  the  computer 


program  operates.  This  example  Is  simplified  In  the  sense  that  only 
two  points  were  used  to  evaluate  the  numerical  analogue  of  the 
emission  Integral.  The  actual  computer  program  will  use  more  points 
A  listing  of  the  computer  program  Is  presented  In  Appendix  A.  The 
program  Is  well  commented  with  Intent  of  being  self-documenting. 


SECTION  IV 


USERS  GUIDE 


Sections  I,  II,  and  III  provide  the  theoretical  basis  for  the 
computer  program  LARTAP.  The  objective  of  this  section  Is  to  provide 
the  user  with  a  guide  of  how  to  use  the  current  version  of  LARTAP. 

The  word,  current.  Is  significant  because  one  of  the  major  guidelines 
utilized  during  the  development  of  LARTAP  was  that  the  program  be 
simple.  The  objective  of  this  simplicity  was  to  allow  for  easy 
modification.  Thus,  the  program  may  be  tailored  with  little  effort  to 
accommodate  any  project  that  requires  the  analysis  of  low  altitude 
radiation  transmission. 

In  order  to  use  LARTAP  the  user  must  place  the  required  program 
Inputs  In  a  disk  file.  This  Is  accomplished  by  utilizing  the  editor 
on  the  user's  computer.  The  disk  file  must  be  named  DATA.  On  each 
line  of  the  file  the  user  must  specify  the  numerical  data  as  presented 
In  the  following  paragraphs.  The  program  utilizes  free  format  read 
statements.  Thus  each  numerical  value  must  be  followed  by  a  comma 
except  for  the  last  piece  of  data  on  a  line. 

On  line  1  of  the  file  DATA,  3  values  must  be  specified,  as 
foil ows : 

MODEL,  GRASL ,  NPT 3 

MODEL  -  An  Integer  that  determines  the  default  atmospheric 

thermodynamic  data  to  be  used  If  the  user  desires.  The  uses 
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of  MODEL  w I  I  I  become  apparent  when  one  enters  the  data  for 
line  3.  The  possible  value  of  MODEL  are: 

MODEL  Default  Atmospheric  Model 

1  Tropical 

2  Midlatitude  Summer 

3  Midlatitude  Winter 

4  Subarctic  Summer 

5  Subarctic  Winter 

6  1962  U.S.  Standard 

Further  documentation  of  these  atmospheric  models  can  be  found 
In  reference  2. 

GRASL  -  A  real  number  that  specifies  the  local  ground  altitude 
measured  relative  to  sea  level  In  kilometers. 

NPTS  -  An  Integer  that  specifies  the  number  of  points  to  be  entered 

after,  and  Including,  line  3  to  define  the  thermodynamic  state 
of  the  atmosphere  as  a  function  of  altitude.  Note  the  larger 
NPTS,  the  more  accurate  the  results  given  by  LARTAP.  A  good 
value  for  NPTS  would  be  10. 

On  line  2  of  the  file  DATA,  3  values  must  be  specified,  as 
foil ows : 

HI,  SMAX,  ANGLE 

HI  -  A  real  number  that  specifies  the  altitude  of  the  radiative 


source  measured  relative  to  the  local  ground  altitude  (GRASL) 
In  k l  I  ometers .. 

SMAX  -  A  real  number  that  specifies  the  distance  from  the  source  to 
the  observer  along  a  straight  line  between  the  source  and  the 
observer.  In  kilometers. 

ANGLE  -  A  real  number  that  specifies  the  angle  between  a  vertical  line 
and  the  line  between  the  source  and  the  observer.  In  degrees. 
This  variable  Is  represented  In  Figure  25  of  Section  III  by 
the  angle  &  . 

On  line  3  through  line  3+NPTS  the  user  specifies  the 

thermodynamic  state  of  the  atmosphere  as  a  function  of  altitude.  On 

v 

each  of  these  lines,  the  user  must  specify  the  following  data: 

ZGR,  P,  T,  RH,  WO,  VIS,  I  HA 

ZGR  -  A  real  number  that  specifies  the  altitude.  In  kilometers, 

relative  to  ground  level  at  which  the  rest  of  the  data  on  the 
I  I ne  I s  meas ured . 

P  -  A  real  number  that  specifies  the  total  atmospheric  pressure. 

In  millibars,  at  altitude  ZGR.  A  value  of  -1000.  will  cause  a 
use  of  default  data  as  specified  by  the  value  of  MODEL. 

T  -  A  real  number  that  specifies  the  atmospheric  temperature.  In 
degrees  C,  at  altitude  ZGR.  A  value  of  -1000.  will  cause  a 
use  of  default  data  as  specified  by  the  value  of  MODEL. 

RH  -  A  real  number  that  specifies  the  relative  humidity  of  the 

atmosphere.  In  percent,  at  altitude  of  ZGR.  A  value  of  -1000. 


will  cause  a  use  of  default  data  as  specified  by  the  value  of 
MODEL. 

WO  -A  real  number  that  specifies  the  density  of  ozone  In  the 

atmosphere, I n  grams  per  cubic  meter,  at  altitude  ZGR.  A  value 
of  -1000.  will  cause  a  use  of  default  data  as  specified  by  the 
value  of  MODEL. 

VIS  -  A  real  number  that  specifies  the  value  of  the  metero I og I ca I 

range.  In  kilometers,  at  altitude  ZGR.  A  value  of  -1000.  will 
cause  a  use  of  default  value  as  specified  by  the  value  of  IHA. 
IHA  -  An  Integer  that  specifies  the  type  of.  particulate  In  the 

atmosphere  region  between  ZGR  and  the  next  specified  altitude 
above  ZGR.  The  possible  values  of  IHA  are  as  follows: 

IHA  Type 

0  No  particulates  present 

1  Rural  with  a  default  metero log  I  cal  range  of  23  km 

2  Rural  with  a  default  metero I og I ca I  range  of  5  km 

3  Maritime  with  a  default  metero I og I ca I  range  of  23  km 

4  Maritime  with  a  default  metero I og  I  ca I  range  of  5  km 

5  Urban  with  a  default  metero I og l ca I  range  of  5  km 

8  Advectlon  Fog  with  a  default  metero I og I ca I  range  of  0.5  km 

9  Radiation  Fog  with  a  default  metero I og I ca I  range  of  0.5  km 

The  Rural  type  Is  Intended  to  represent  the  particulate 
conditions  found  In  continental  areas  which  are  not  directly 
Influenced  by  urban  particulate  sources. 

The  Maritime  type  applies  to  oceanic  regions  and  coast  line 


areas.  The  particulates  result  from  sea-salt  part  »,  ,■ 
produced  by  the  evaporation  of  sea-spray  droplets. 

The  Urban  type  applies  to  Industrial  and  urban  areas  •  •«  -  • 
majority  of  particulates  result  from  combustion  processes. 

When  the  atmosphere  becomes  nearly  saturated  with  water  vapor, 
fog  can  form.  Saturation  of  the  air  can  occur  as  a  result  of  two 
different  processes;  the  mixing  of  air  masses  with  different 
temperatures  (advectlon  fog)  or  by  cooling  of  the  air  to  the  point 
where  Its  temperature  approaches  the  dew-point  (radiation  fog )  . 

For  thick  fogs  where  V  Is  less  than  200m  the  extinction 
coefficient  Is  virtually  Independent  of  wave  length.  These  conditions 
are  modeled  best  by  the  advectlon  fog.  For  moderate  fogs, 

200m  V  1000m,  the  radiation  fog  should  be  used.  For  thin  fogs, 

1km  V  2km,  a  maritime,  urban  or  rural  particulate  model  should  be 
used  with  a  99$  relative  humidity. 

It  should  be  mentioned  that  a  cloud  Is  nothing  more  than  a  thick 
fog  that  does  not  contact  the  ground. 

Since  the  density  of  ozone  at  low  altitudes  Is  quite  low  It  Is 
suggested  that  no  special  effort  be  undertaken  to  measure  ozone 
density.  However,  for  completeness.  It  Is  recommended  that  ozone  be 
Included  utilizing  one  of  the  program  default  model  capabilities. 

In  order  for  LARTAP  to  operate  correctly,  NPTS  must  have  a  value 
of  at  least  two.  This  means  that  the  thermodynamic  state  of  the 
atmosphere  must  be  specified  at  a  minimum  of  two  altitudes,  (two 
values  of  ZGR).  These  two  altitudes  have  to  bound  the  atmospheric 
region  of  Interest.  For  example  ZGR  =  0.0  and  ZGR  =  0.7  where  0.0 
represents  the  minimum  altitude  of  Interest  relative  to  ground  level 
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•’  represents  the  maximum  altitude  of  Interest  relative  to  ground 

o  i  S  . 

The  output  from  the  current  version  of  LARTAP  Is  In  the  form  of 
three  disk  files.  The  names  of  the  disk  files  are  WAVN,  TRAS,  and 
RADA.  The  file  WAVN  contains  the  wave  numbers  for  the  spectral  region 
studied  by  LARTAP.  WAVN  contains  one  wave  number  on  each  line  of  the 

file.  The  file  TRAS  contains  the  transmittance  for  the  atmospheric 
configuration  specified  In  DATA.  TRAS  contains  one  value  of 
transmittance  on  each  line.  The  value  of  transmittance  on  a 
particular  line  of  TRAS  corresponds  to  the  wave  number  on  the  same 
line  In  WAVN.  The  file  RADA  contains  the  Intensity  of  radiation 
emitted  by  the  atmosphere  between  the  source  and  the  observer  for  the 
atmospheric  configuration  specified  In  DATA.  The  value  of  emitted 
Intensity  on  a  particular  line  In  RADA  corresponds  to  the  wave  number 
on  the  same  line  In  WAVN. 

The  spectral  region  studied  by  the  current  version  of  LARTAP  Is 
from  a  wave  number  of  350  cm  to  a  wave  number  of  40,000  cm  In 
Increments  of  10  cm 

In  table  1,  five  sample  Input  data  flies  are  presented.  Each  of 
the  five  sample  data  files  shown  In  table  I  Is  preceded  by  a  title 
Indicating  the  values  of  HI,  SMAX,  and  ANGLE  for  that  particular  data 
file.  These  titles  are  not  part  of  the  data  files.  LARTAP  output  for 
these  five  sample  data  flies  Is  presented  graphically  In  figure  26 
through  figure  40. 

The  program  also  requires  seven  other  Input  data  flies.  These 
files  are  not  controlled  by  the  user.  They  supply  the  necessary 
spectral  data  that  was  discussed  In  section  III.  The  reasons  for 
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storing  this  spectral  data  In  flies  Is  as  follows. 

In  order  to  facilitate  this  possible  conversion  of  this  program 
to  a  micro-computer  for  use  in  the  field.  It  was  determined  that  the 
large  spectral  data  based  would  have  to  be  stored  In  peripheral 
equipment.  Thus  these  seven  spectral  data  files  were  Incorporated  to 
emulate  such  peripheral  equipment. 


TABLE  1 


SAMPLE  DATA  FILES 


M1  =  l  O  KM  SMAX^l.O  KM  ANGLE--0  DEGREES 


1000.  0/  0 
1000.  0,  0 
1000.  o,  o 
1000.  o,  o 
1000.  o.  0 


Hl-1.0  KM  SMAX-1.  0  KM  ANGLE=60  DFGRI-.E3 

1-  0.  000.  5 
1  0,  1.0,  60.  000 

0  000,  -1000.0,  -1000.0,  -1000.0,  -1000  0, 

0  500,  -1000  0,  -1000.0,  -1000.0,  -1000  o, 

1  000,  -1000  0.  -1000.0,  -iOOO.  0,  -1000.0, 

1  500.  -1000  0,  -1Q00.0,  "1000.0,  -IOOO  o, 

2  000,  -1000  O,  -1000.0,  -1000.0,  -1000  o, 


HI  *  1 .  O  KM  SMAX=»  1  0  KM  ANGLE-90  DEGREES 


1, 0  000,  5 
1  0,  1  0,  90  000 


0  000. 

-1000.  o. 

-1000.  0, 

1000.  0, 

-1000.  0, 

1000.  0, 

0.  500> 

-1000  o, 

-1000  0. 

-1000  o. 

-1000.  o, 

-1000  0, 

1  000/ 

-1000.  0, 

-1000  0, 

'1000.  o, 

-1000  o. 

-1000.  0, 

1 .  500  < 

-1000  o, 

-1000.  0, 

-1000.  0, 

1000  0, 

-IOOO.  0, 

2.  000 j 

-1000  0, 

-1000.  0, 

-1000.  0, 

-1000  0, 

1000.  0, 

HI  =  ! .  0  KM  3MAX=1.0  KM  ANGLE— 1 20  DECREES 


1000  0,  0 
1000.  O,  0 
1000.0,  G 
1000  0,  0 
1000  0,  0 


SI 


1. 0.  000,  5 

1  o,  1  o.  120.  000 

0.000,  -1000.0,  -1000.0,  "1000.0,  -1000  o,  - 

0.  500,  -1000  0,  -1000.  0,  -1000.  0,  -1000  0, 

1.000,  1000.0,  -1000.0,  -1000.0,  -1000  o,  - 

1.500,  -1000.0,  -1000.0,  "1000.0,  -1000  0,  • 

2.000,  -1000.0,  -1000.0,  -1000  0,  -1000.0, 


-1000.  O,  0 
-IOOO  O,  0 
-1000.  O,  0 
-IOOO  O,  o 
-1C00.  O,  u 


1,  0.  000,  5 

1  o,  1.  0,  0.  000 

0.  000,  -1000.  0,  -1000.  0,  -1000.  0,  -1000.  0,  - 

0.500,  -1000.0,  -1000.0,  -1000.0,  -1000.0,  - 

1.  000,  -1000.  0,  -1000.  0,  -1000.  0,  -1000.  0,  - 

1.500,  -1000.0,  -1000.0,  -1000.0,  -1000.0,  - 

2.  000,  -1000  0,  -1000.  0,  -1000.  O,  -1000.  o, 


TABLE  1  (CONTINUED) 


Hl*l.  0  KM  SMAX- 1 .  O  KM  ANGLE" 1 30  DEGREES 


1.  0  000,  5 
1.0.1.  0,  120.  000 


0.  000, 

-1000.  0. 

-1000.  0, 

-1000.  0. 

-1000.  0. 

-1000.  0, 

0 

0  500, 

-1000.  0, 

-1000.  0, 

-1000.  0, 

-1000. 0, 

-1000.  0, 

0 

1.  000, 

-1000.  0, 

-1000.  0. 

-1000.  0, 

-1000.  0, 

-1000.  0. 

0 

1.  500, 

-1000.  0, 

-1000.  0, 

-1000.  0, 

-1000.  0, 

-1000.  0, 

0 

2.  000, 

-1000.  0, 

-1000.  0, 

-1000.  0, 

-1000.  0, 

-1000.  0, 

0 

5 


s  vs  *. 


FIGURE  S8 

Hi-1  KM  RANGE-1  KM  ANGLE-0  DEGREES 


FIGURE  31 

Hl-1  KM  RANGE-1  KM  ANGLE-60  DEGREES 


FIGURE  32 

Hl-1  <H  RANGE- 1  KM  ANGLE -90  DEGREES 


UAUENUMBER  -  CM-1  THOUSANDS 


FIGURE  33 

Hl-1  ICM  RANGE- 1  KM  ANGLE-90  DEGREES 
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SECTION  V 


SUMMARY  OF  RADIATION  INSTRUMENTS 

It  can  be  concluded  from  section  IV  that  In  order  to  support  the 
current  version  of  LARTAP,  It  Is  necessary  to  make  four  measurements. 
These  four  measurements  are  pressure,  temperature,  relative  humidity 
and  metero I og I ca I  range.  As  was  discussed  In  section  IV,  since  the 
density  of  ozone  at  low  altitudes  Is  small.  It  Is  suggested  that  no 
special  effort  be  made  to  measure  ozone  density.  However,  for 
completeness.  It  is  recommended  that  ozone  be  Included  utilizing  one 
of  the  program  default  model  capabilities. 

The  pressure  measurement  can  be  made  utilizing  a  barometer.  The 
temperature  and  relative  humidity  measurements  can  be  made  utilizing 
an  AN/TMQ-11  Temperature  Dew  Point,  finally,  the  metero I og I ca I  range 
measurement  can  be  made  utilizing  an  AN/GMQ-10  TransmI ssometer .  This 
Instrumentation  Is  available  at  most  AFB  Installations. 

The  atmospheric  components  that  effect  radiation  transmission  are 
H  40 ,  0JL,  CO,,  CH4,  N*0,  CO,  0,,  Na  and  particulates.  The  question  now 
arises  as  to  how  each  of  these  components  Is  accounted  for  by  the  four 
required  measurements.  The  H40  component  Is  accounted  for  by  the 
relative  humidity  measurement.  The  03  component  would  ordinarily  be 
accounted  for  by  direct  measurement,  however,  as  was  Just  rilscussed.lt 
Is  recommended  that  default  data  be  utilized  so  no  measurement  Is 
necessary.  The  other  gaseous  components  and  the  particulate  component 
are  accounted  for  by  the  use  of  two  assumptions.  First,  In  order  to 
account  for  the  remaining  gaseous  components  It  Is  assumed  that  the 
mole  fraction  concentration  of  each  of  these  components  Is  a  world 
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wide  constant.  Secondly,  as  was  discussed  In  section  III,  by  assuming 
standard  size  distributions,  the  particulate  component  can  be 
accounted  for  by  utilizing  the  metero I og I ca I  range.  These  two 
assumptions  are  Inherent  In  the  data  base  obtained  from  LOWTRAN  for 
use  In  LARTAP.  According  to  reference  2  these  assumptions  are  valid 
for  the  Intended  uses  of  LOWTRAN.  However.,  LARTAP  will  be  used  by 
AFFTC  at  Edwards  AFB  where,  state  of  the  art,  radiation  systems  are 
being  tested.  Thus  the  question  arises  as  to  the  accuracy  of  these 
assumptions  for  the  Intended  uses  of  LARTAP.  This  question  can  only 
be  answered  by  comparing  results  predicted  by  LARTAP  with  these 
obtained  from  test  data.  Should  disagreement  arise  then  the  LARTAP 
data  base  would  have  to  be  expanded  through  detailed  testing  and  more 
detailed  measurements  would  have  to  be  taken  to  support  LARTAP. 

CSUN  has  surveyed  the  equipment  available  at  Edwards  AFB  to 
determine  the  currently  existing  capability  to  support  such  a  detailed 
testing  program.  The  present  Inventory  at  Edwards  AFB  contained 
Instrumentation  to  measure  HA0,  05,  and  CO  concentrations.  The 
following  Is  an  outline  of  this  Instrumentation. 

TMQ-r  1 1  : 

This  Is  standard  equipment  at  most  AFB  Installations  capable  of 
Indicating  relative  or  absolute  water  content  of  the  air. 

Mooli, or  Labs  8410^.  Qzono-AnaJ yaar 

The  concentration  of  ozone  Is  measured  by  a  chemiluminescent  detector. 
The  ozone  sample  and  ethylene  are  mixed  In  a  reaction  chamber.  The 
ozone  reacts  with  the  ethylene  to  form  activated  aldehyde.  The 
aldehyde  then  emits  Infrared  with  an  Intensity  directly  proportional 
to  the  ozone  concentration. 
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Full  Scale  Ranges  -  0.1,  0.5,  1.0,  5.0,  10.0  ppm. 

Precision  -  +  0.339$  full-scale 

Minimum  Sensitivity  -  2  ppb. 

Monitor  Labs  8310  -  Carbon  Monoxide  Analyzer*. 

Ambient  air  Is  pumped  In  and  split  Into  two  paths.  One  path  flows 
through  a  sample  cell  where  an  Infrared  beam’s  Intensity  Is  measured. 
The  higher  the  CO  content,  the  greater  the  attenuation  of  the  IR  beam 
The  second  flow  path  travels  through  a  catalytic  converter,  which 
changes  all  the  CO  to  C0a ,  and  then  flows  to  a  second  sample  cell. 
This  second  cell  Is  used  as  a  reference,  and  the  difference  In  the  IR 
attenuations  Is  used  to  display  the  CO  concentration. 

Full  Scale  Ranges  -  50,  100  ppm. 

Precision  -  +  0.1  ppm. 

Minimum  Sensitivity  -  0.1  ppm. 

Edwards  AFB  also  has  a  good  support  Instrument  available  In  the 
following: 

Han  I  tflr_-1.flfas-93QQj _ Data  Logger 

This  Instrument  Is  a  general  purpose,  microprocessor  controlled,  data 
logger.  Analog  Inputs  are  attached  to  any  of  the  40  channels,  where 
the  signal  Is  digitized  and  stored.  The  data  flow  Is  programmed  via 
the  fron^  key  board.  Standard  features  Include:  English  language 
programming  with  alpha-numeric  display,  auto-ranglng  on  Input 
voltages,  thermocouple  printout,  and  alarm  conditions.  Optimal 
features  available  are:  alarm  relay  panel,  external  recording  devices 
special  order  conversion  fables,  and  programmable  alarms. 

No.  of  Channels  -  40 

Channel  Voltage  Ranges  -  10.000V,  3.000V,  300.00mV,  30.000mV. 
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Input  Impedance  -  100  Meg  ohm 
Channel  Scan  Rate  -  1  sec.  to  999  hrs. 

Printer  Scan  Rate  -  1  sec.  to  999  hrs. 

Optional  External  Recording  Devices  - 
Tel etype 

Serial  Terminals  and  Printers 
Paper  Tape  Punches 
Incremental  Mag.  Tape 

Synchronous  Mag.  Tape  with  read  after  write  data  check 
Alarms  -  Channel  Overrange  (exceeds  full  scale) 

Thermocouple  burned  out  (open  channel) 

In  an  effort  to  complete  the  required  capability.  Instrument 
manufacturers  were  contacted.  CSUN  was  able  to  acquire  Information 
concerning  Instruments  capable  of  measuring  C0^  concentrations,  CH^. 
concentrations  and  particulate  concentrations  (and  size  distribution) 
The  following  Is  a  summary  of  this  Instrumentation. 

Monitor.  Labs  33QQ/3400i  Carbon  Monox  I de/D I  ox  I de  Meter 
This  Instrument  utilizes  a  single  beam  Infrared  source  In  conjunction 
with  a  mechanical  chopper.  The  3300  Is  single  component  and  3400  dual 
component . 

Full  Scale  Range  C04  -  400  ppm. 

Full  Scale  Range  Co  -  1000  ppm. 

Precision  -  +1$ 


Price  3300/3400  -  $2,975/$4,500 


Logger.  The  Model  400  measures  hydrocarbons  utilizing  flame 
ionization  detection.  The  sample  Is  Injected  Into  a  burning 
hydrocarbon  flame  where  the  hydrocarbons  become  Ionized.  A  voltage  Is 
applied  between  the  burner  Jet  and  a  collection  plate  that  results  In 
Ion  migration.  The  Induced  current  Is  proportional  to  the  hydrocarbon 
concentrat I  on . 

*F u I  I  Scale  Range  -  1,  10,  100,  1000  ppm. 

Precision  -  ±\%  of  full  scale 

Output  -  lOmV,  lOOmV,  IV,  5V  (selectable) 

Price  -  $4,600 
Recommended  Accessories: 


Part  Ho. 

OftSET  IptJfflO. 

Pr.lfifi 

630779 

Air  Pump 

$  410 

630157 

Auto  Fuel  Shut-off 

$  145 

Sollnold  for  Shut-off 

$  57 

*WIth  pure  hydrogen  fuel  accessory 

It  Is  possible  to  obtain  more  detailed  measurements  of 
hydrocarbons  by  using  a  gas  chromatograph.  The  chromatograph  can 
break  down  the  different  hydrocarbon  component,  but  such  an  Instrument 
must  be  designed  for  the  particular  tests,  and  chromatographs  are  in 
the  $20,000  price  range. 

Various  particle-sizing  Instruments  were  evaluated  for  measuring 
particulate  concentration  and  size  d I str I  but  I  on .  Those  capable  of 
measuring  particles  of  size  less  than  10  micrometers  with  reasonable 
accuracy  and  ease  ot  operation  were  selected.  Those  Instruments  found 
to  satisfy  these  requirements  are  presented  below. 


Anderson.  Samp  I ers  Inc.; _ Low-Pressure  Impactor 

The  Anderson  LPI  utilizes  the  well-known  Inertial  Impactor  principle 
of  the  particles  In  the  gas  flow;  the  higher  Inertia  of  the  heavier 
phase  can  cause  the  particle  trajectories  to  differ  from  the  gas  flow 
streamlines.  By  staging  plates  with  nozzles  of  decreasing  diameter, 
size  differentiation  Is  obtained;  as  the  gas  velocity  Increases, 
successively  smaller  particles  will  Impact  on  the  dish  beneath  the 
nozzle.  The  collected  particles  are  transported  to  the  laboratory  for 
ana  lysis. 

The  price  of  the  complete  system  Is  $4,875.  The  disadvantage  Is 
accuracy  may  be  lost  unless  extreme  care  Is  exercised  In  handling  and 
transporting  the  collected  particles. 

Cal.J..for.nla..Msasur9m9nts . Ins. ; _ As  loss.]...  Part  Isle-Anal  y.z  er_Mg.de] 

EQr-Z 

The  PC-2  Is  a  10  stage  cascade  Impactor  with  quartz  crystal 
microbalance  mass  monitors  In  each  stage.  It  utilizes  the  Inertial 
Impactor  principle  as  In  the  Anderson  L.P.I.,  however,  the  PC-2 
processes  a  list  of  special  features  which  make  It  most  attractive  for 
Its  accuracy  and  ease  of  operation.  The  most  prominent  of  Its 
features  Is  the  employment  of  highly  sensitive  piezo  electric  crystal 
mass  monitor  In  each  stage  which  provides  real  time  mass  concentration 
I  n  f  or/nat  I  on .  The  adhesive  coated  surface  of  the  crystals  retain  the 
samples  for  further  evaluation. 

Other  features  Include  built-in  crystal  saturation  Indicator  and 
bullt-tn  sample  time  stopwatch. 

The  estimated  price  for  the  complete  system  $15,000. 


CONCLUSION  - 


As  Infrared  radiation  technology  becomes  more  and  more  a  part  of 
our  defense  capabilities,  the  need  arises  to  be  able  to  understand  and 
predict  the  effects  of  the  earth's  atmosphere  on  the  transmission  of 
Infrared  radiation.  The  need  Is  especially  great  for  the  low  altitude 
regions  of  the  earth's  atmosphere. 

LOWTRAN  offers  one  possible  solution  to  this  problem.  However,  as 
LOWTRAN  has  evolved  It  has  become  a  very  complicated  computer  program. 
Thus  a  new  program,  LARTAP,  has  been  developed.  LARTAP  utilizes  the 
large  empirical  data  base  from  LOWTRAN,  however,  this  new  program  Is 
geared  specifically  towards  low  altitude  radiation  problems,  thus  the 
resulting  computer  code  Is  much  simpler  than  LOWTRAN. 

The  net  result  Is  that  LARTAP  offers  a  good  staitlng  point  for  the 
development  of  any  computer  program  to  support  a  project  that  requires 
Information  concerning  low  altitude  radiation  transmission. 

One  of  the  specific  goals  of  this  project  was  to  attempt  to  make 
LARTAP  a  tool  that  could  be  used  In  the  field  to  provide  analysis 
capability  at  a  test  site.  This  would  require  the  conversion  of  the 
program  to  a  m 1 cro-computer .  In  order  to  accomplish  this  conversion, 
the  large  empirical  data  base  utilized  by  LARTAP  would  have  to  be 
stored  on  peripheral  equipment.  Thus  In  the  current  version  of  LARTAP 
this  data  base  has  been  stored  In  disk  flies  to  emulate  such 
peripheral  equipment. 

A  short  coming  of  both  LOWTRAN  and  LARTAP  Is  the  neglect  of 
scattering  Into  the  line  of  sight  of  radiation  from  sources  other  than 
the  one  of  Interest.  It  Is  felt  that  this  could  be  quite  significant 
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for  the  scattering  of  solar  radiation.  It  Is  suggested  that  an  effort 
should  be  under  taken  to  model  this  effect.  The  resulting  model  could 
easily  be  Incorporated  Into  LARTAP. 

It  Is  suggested  that  the  LARTAP  program  be  tested.  If  the  program 
should  not  agree  with  test  data  then  detailed  testing  would  have  to  be 
performed  to  modify  the  data  base  of  LARTAP. 

Although  the  primary  concern  of  this  project  has  been  the  modeling 
of  radiation  transmission,  CSUN  has  also  considered  the  possibility  of 
direct  measurement  of  radiation  transmission  at  a  test  site. 

More  direct  techniques  of  measuring  radiation  transmission  can  be 
accomplished  by  using  a  radiometer.  A  wide  range  radiometer,  such  as 
the  Molectron  PR200,  (see  table  2),  Is  sensitive  from  the  ultraviolet 
to  the  far  Infrared  ranges.  By  using  Interchangeable  filters,  the 
proper  frequency  range  can  be  obtained. 

For  example,  the  radiometer  Is  mounted  on  a  GMD-2  Balloon  Tracker 
to  track  a  radiation  source  that  Is  rising  through  the  atmosphere  on  a 
balloon.  The  light  source  has  to  be  stable  and  also  flash  on  and  off 
so  that  changing  background  radiation  can  be  measured. 

The  measurements  from  this  test  could  be  combined  with  the  theory 
of  sections  I  and  II  of  this  report  utilizing  a  numerical  filtering 
technique.  The  result  would  be  an  accurate  model  of  radiation 
transmission  at  the  test  site. 

In  conclusion,  a  simple  computer  program,  LARTAP,  has  been 
developed.  This  program  needs  to  be  tested.  The  simplicity  of  the 
program  makes  It  essy  to  modify  to  suit  users  requirements.  Should 
detailed  testing  be  required  to  Improve  program  accuracy,  currently 
available  Instrumentation  has  been  evaluated.  It  has  also  been 


TABLE  2 


Mo  I ectron.  PR2QQ  Pyroelectric  Radiometer 

Detector  Type  -  Circular  lithium  tantalate  pyroelectric 

Spectral  Range  -  0.2  to  20  +3f 

0.2  to  40  <u.m  ±51 

40  to  500  «.m  uncallbrated 

Ranges  (Full  Scale)  -  2  ^W,  20  200  *W,  2  mW,  20  mW, 

System  Response  Time  -  Full  Scale  Range  Fast  Response 

10  sec. 

3  sec. 

2  sec. 

1  sec. 


2  uW 
20  icW 
200  ziW 
2,20,200  mW 


crystal 


200  mW/cm 
Slow  Respon 
100  sec. 
100  sec. 
100  sec. 
100  sec. 


Absolute  Maximum  Irradlance  -  50  Watts/cm 
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APPENDIX  A 


Listing  of  Computer 
Program 
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5 


1 


LOW  ALTITUDE  RADIATION  TRANSMISSION  ANALYSER  PROGRAM 

L.  A.  R.  T.  A.  P. 

MAIN  DRIVER  PROGRAM 

COMMON  /BLK1 /  ZGR  ( 20 ) ,  P ( 20 ) , T ( 20 ) ,  RH ( 20 ) , WH<  20 ) »  WO ( 20 ) .  VIS ( 20  > 
COMMON  /BLK2/  HI  ,  ANGLE  ,  SMAX 

COMMON  /BLK3/  IHA(20),  NPTS 

COMMON  /SPEC1/  Cl ARR  <  2580 ) ,  C2ARR ( 1575).  C3ARR  <  540 ) ,  C4ARR ( 133 ) 
COMMON  /SPEC2/  C7ARR ( 280,  5 > »  C8ARR(102.2)  ,C9ARR(131) 

DIMENSION  Wl(10),  W2(10),  REL(20),  WLAY(20, 10),  CVE(IO),  CVA(IO) 
DIMENSION  EXTN ( 20,  10),  ABSP(20,  10),  EMIS(20) 

DIMENSION  TEMP (20),  TRANE (20),  TRANA(20) 

DIMENSION  IHAZ ( 20 ) 

OPEN  OUTPUT  FILES 
0PEN(UNIT=6,  FILE='WAVN' ) 

0PEN(UNIT=7,  FILE=?'TRAS ' > 

OPEN ( UNIT=8,  FILE* 'RADA  ' ) 

READ  SPECTRAL  DATAFILES 
OPEN ( UNIT=31 ,  FILE™ 'C 1FILE ' ) 

0PEN(UNIT=32,  FILE* 'C2FILE ' ) 

OPEN ( UNIT=33,  FILE* 'C3FILE '  > 

OPEN < UN I T *34,  FILE* 'C4FILE ' ) 

OPEN ( UN I T =37 ,  FILE* 'C7FILE ' ) 

0PEN<UNIT=38,  FILE* 'C8FILE ' ) 

0PEN(UNIT=40,  FILE* 'C9FILE ' ) 

H20  SPECTRAL  DATA 
DO  31  1  =1,2580 
READ (31,  *)  CIARR(I) 

CONTINUE 

UNIFORM  GASES  SPECTRAL  DATA 
DO  32  1  =  1,  1575 
READ (32,*)  C2ARR ( I ) 

CONTINUE 

OZONE  I.R.  SPECTRAL  DATA 
DO  33  1  =  1,  540 
READ (33,*)  C3ARR ( I ) 

CONTINUE 

N2  SPECTRAL  DATA 
DO  34  1  =  1,  133 
READ (34,*)  C4ARR ( I ) 

CONTINUE 

AEROSOL  SPECTRAL  DATA 
DO  37  1=1, 280 

READ  (  37,  *  )  C7ARR  (  1 ,  1  ) ,  C7ARR  (1,2).,  C7ARR  (1,3),  C7ARR  (1,4),  C7ARR  (1,5) 
CONTINUE 

OZONE  U.  V.  SPECTRAL  DATA 
DO  38  1  =  1,  102 

READ  (38,*)  C8ARR  ( 1 ,  1  )  ,  C8ARR  (1,2) 

CONTINUE 

H20  CONTINUUM  (4  UM)  SPECTRAL  DATA 
DO  39  1  =  1,  131 
READ  (40,*)  C9ARR  ( I ) 

CONTINUE 


DATA  INITIALIZE 
CALL  PROPUT 

CALCULATE  COMPONENT  AMOUNT  INTEGRALS 

LY=0 

S1=0.  OOO 

-  CALCULATE  COMPONENT  AMOUNTS  AT  START  OF  PATH 
CALL  EQUABS(S1,  IHAZE1, Ml,  REL1.  TEM1 ) 

INTEGRAL  EVALUATION  LOOP 

100  LY=LY+1 

S2=SMAX  +  1.  0E-04 

CHECK  FOR  HORIZONTAL  PATH 

IF  (ANGLE  .  GT.  89.3  .AND.  ANGLE  .  LT.  90.5)  GO  TO  120 
FIND  CLOSEST  DATA  POINT  -  DEFAULT  IS  SMAX 
DO  110  1  =  1,  NPTS 
DSMIN=S2— SI 

TEST* ( ( ZGR ( I )— HI ) /COS ( ANGLE*3.  14159/180.  O) )-Sl 
IF  (TEST  . GT.  0.  0001  .  AND.  TEST  .  LT.  DSMIN)  S2=TEST+S1 
110  CONTINUE 

C  CALCULATE  COMPONENT  AMOUNT  INTEGRALS  FOR  LAYER  LY 

120  CALL  EQUABS(S2<  IHAZE2.  W2, REL2<  TEM2) 

DO  130  1*1.9 

WLAY(LY.  I )  =  (W2( I )+Wl ( I ) >*(S2-S1 )/2.  O 

IF  (W1(I)  .  EQ.  0.000  .OR.  W2(I>  .  EQ.  0.000)  GO  TO  130 
AUG=W2(I)/W1(I) 

WLAY(LY. I ) =W 1 ( I ) * ( S2— S 1 ) 

IF  (AUG  .  GT.  0.99  .AND.  AUG  .  LT.  1.01)  GO  TO  130 
WLAY(LYj I ) = ( W2 ( I ) -W 1 ( I ) ) * ( S2— S 1 ) / ALOG ( AUG ) 

130  CONTINUE 

C  CALCULATE  AVERAGE  RELATIVE  HUMIDITY  FOR  LAYER  LY 

REL(LY) =(REL1+REL2) /2.  0 

C  DETERMINE  AEROSOL  TYPE  FOR  LAYER  LY 

IF  (ANGLE  .  LE.  90.0)  IHAZ(LY)=IHAZE1 
IF  (ANGLE  .  GT.  90.0)  IHAZ(LY)  =  IHAZE2 
C  STORE  TEMPERATURE  FOR  EMISSION  CALCULATION 

TEMP ( LY ) =TEM 1 

C  INCREMENT  DATA 

S1=S2 
REL1=REL2 
IHAZE1=IHAZE2 
TEM1=TEM2 
DO  250  1*1,  9 
W1 ( I )=W2( I ) 

250  CONTINUE 

C  CHECK  FOR  END  OF  PATH 

IF  ( S2  .  LT.  SMAX)  GO  TO  100 

LYT0T=LY+1 

C  STORE  TEMPERATURE  AT  SMAX 

TEMP ( LYTOT ) =TEh2 
C  SPECTRAL  LOOP 

DO  400  IV=350,  40000,  10 

C  PATH  INDEPENDENT  SPECTRAL  COEFFICIENTS 

CALL  C1DTA(CVE( 1 ) ,  IV) 

CVA ( 1 ) =CVE ( 1 ) 

CALL  C2DTA ( CVE ( 2 ) ,  IV) 
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^  108 

CVA(2)“CVE(2) 

►V  109 

CALL  C3DTA( CVE ( 3 ) <  IV) 

110 

CVA  <  3 ) “CVE  <  3 ) 

L  ■  111 

CALL  C4DTA ( CVE ( 4 ) , IV) 

v  112 

CVA ( 4 ) “CVE ( 4 ) 

CALL  C5DTA(CVE<  5) >  IV) 

"  114 

CVA( 5 ) =CVE(  5 ) 

115 

CALL  C6DTA<CVE< 6 ) ,  IV) 

•  116 

CVA<6)=0.  000 

■  ii7 

CALL  £8DTA(CVE(8),  IV) 

•V  118 

CVA ( 8 ) =C VE  <  8 ) 

|j|  1 19 

CALL  C9DTA(CVE(9), IV) 

P  '  120 

CVA ( 9 ) “CVE ( 9 ) 

'  121 

c 

ZERO  EXTINCTION  AND  ADSORPTION  ARRAYS 

122 

DO  265  J“l,  9 

123 

EXTN(LYTOT. J)=0.  000 

124 

ABSP (LYTGT,  J)=0.  000 

125 

265 

CONTINUE 

■  126 

C 

CALCULATE  TRANSMITTANCES  FROM  EACH  LAYER 

TO  OBSERVER 

127 

C 

NOTE  -  WORK  FROM  OBSERVER  BACK  TO  SOURCE 

128 

I=LYTOT 

129 

DO  300  K=l»  LYTOT 

1 30 

C 

CALCULATE  TRANSMITTANCE  AT  WAVENUMBER  IV 

(EXTINCTION) 

^  131 

C 

H20 

M  132 

TRANE ( I  )  =  1.  00/EXP  (0.  06S9*EXTN(  I,  1  )**0.  555) 

133 

C 

UNIFORM  GASES 

■>  1 34 

TRANE ( I ) =TRANE ( I ) /EXP  < 0.  0689*EXTN( I,  2)**0.  555) 

135 

C 

OZONE  I.  R. 

V.  136 

TRANE ( I ) “TRANE ( I ) /EXP ( 0.  0643*EXTN ( I >  3 ) **0.  666 ) 

£  137 

C 

N2  CONTINUUM 

I  138 

TRANE ( I )=TRANE< I ) /EXP < EXTN ( I, 4) ) 

“  139 

C 

H20  CONTINUUM  (10  UM) 

140 

TRANE  ( I  )  =»TRANE  ( I )  /EXP  (EXTN(  I ,  5)  ) 

141 

C 

MOLECULAR  SCATTERING 

142 

TRANE ( I ) =TR ANE ( I ) /EXP ( EXTN ( 1 ,  6 ) ) 

143 

C 

AEROSOL 

Ml  144 

TRANE ( I ) “TRANE < I ) /EXP  <  EXTN (1,7)) 

145 

C 

OZONE  U.  V. 

■\  1 46 

TRANE ( I ) “TRANE ( I ) /EXP ( EXTN ( 1 ,  8 )  ) 

£  147 

C 

H20  CONTINUUM  (4  UM) 

148 

TRANE ( I ) “TRANE ( I ) /EXP ( EXTN( 1 ,  9) ) 

;>  i  49 

C 

CALCULATE  TRANSMITTANCE  AT  WAVENUMBER  IV 

(ABSORPTION) 

't  150 

C 

H20 

■  151 

TRANA(I)  =  1.  00/EXP (0.  0689*ABSP(I,  1 )**0.  555) 

152 

C 

UNIFORM  GASES 

1 53 

TR ANA  ( I  >  “TRANA  (  I )  /EXP  ( 0.  0689*ABSP  (1,2)  **0.  555 ) 

v-  154 

c 

OZONE  I.  R. 

155 

TRANA ( I ) =TRANA ( I ) /EXP ( 0.  0643*ABSP (1,3) **0.  666 ) 

156 

c 

N2  CONTINUUM 

!|  157 

TRANA ( I ) =TRANA ( I )/EXP (ABSP ( I,  4) ) 

158 

c 

H20  CONTINUUM  (10  UM) 

159 

TRANA ( I ) “TRANA ( I ) /EXP ( ABSP (1,5)) 

V  1 60 

c 

MOLECULAR  SCATTERING 

i6i 

4 

TRANA ( I ) “TRANA ( I > /EXP ( ABSP (1,6)) 

i 

t$ 

82 

V 

%  ,  ■  ^  . 

VI.: : • 1   ..  -A  AhA.-,  i  ■«  i.  ■  ■  , , 

-  - - 

n  n 


C 

C 

C 


C 


275 

300 

C 

325 

350 

C 


375 

C 


400 

1000 


AEROSOL 

TRANA ( I >  =TR ANA ( I ) /EXP ( ABSP ( I » 7) ) 

OZONE  U.  V. 

TRANA( I )=TRANA< I ) /EXP (ABSP ( I,  8) ) 

H20  CONTINUUM  <4  UM) 

TRANA ( I ) = TRANA < I ) /EXP  <  ABSP (1,9)) 

1  =  1-1 

IF  (I  .  EQ.  0)  GO  TO  325 

AEROSOL  SPECTRAL  COEFFICIENTS 
CALL  C7DTA(CVE(7), CVA<7),  IV, REL ( I ) ,  IHAZ  < I ) ) 

CALCULATE  EXTINCTION-ABSORPTION  COFFICIENTS 
FROM  BEGINNING  OF  LAYER  I  TO  OBSERVER 
DO  275  J=l.  9 

EXTN < I ,  J ) =EXTN ( I +1 ,  J ) +CVE ( J ) *WLAY ( I,  J) 

ABSP ( I , J)=ABSP< I  +  l.  J)+CVA( J)*WLAY( I.  J) 

CONTINUE 

CONTINUE 

CALCULATE  EMISSION  ALONG  PATH 
DO  350  1=1 , LYTOT 
EMIS( I )=PLAN(TEMP ( I ) ,  IV) 

CONTINUE 

CALCULATE  RADIANCE  FOR  PATH 
RAD  1=0.  000 
DO  375  1=1, LY 

DRAD=(EMIS< I  +  l )*TRANE< I  +  l ) /TRANA ( I  +  l ) >+<EMIS< I )*TRANE( I ) / TRANA ( I 

RADI=RADI+DRAD*< TRANA <1+1 ) -TRANA ( I ) )/2.  O 

CONTINUE 

PRINT  RESULTS 
WRITE (6, 1000)  IV 
WRITE(7. 1000)  TRANE ( 1 > 

WRITE<8, 1000)  RADI 

CONTINUE 

FORMAT < 1PE12.  5) 

STOP 

END 


83 


C**  PROPUT  *****************  W****************-*-******  ****#■■«••»•*•■«■  *********** 
C 

SUBROUTINE  PROPUT 

C  DATA  PROCESSING  AND  INPUT  ROUTINE 

COMMON  /BLK1  /  ZGR ( 20) , P ( 20 ).T( 20 ) , RH( 20 ) »  WH ( 20 ) , WO  <  20 > , VIS ( 20 ) 
COMMON  /BLK2/  HI  ,  ANGLE  ,  SMAX 

COMMON  /BLK3/  IHA(20),  NPTS 

C  OPEN  INPUT  DATA  FILE 

0PEN<UNIT»5,  FILE-'DATA' > 

C  DEF I ME  CONSTANTS 

RV=4.  6150E-03 

C  PROGRAM  CONTROL  DATA 

READ v  5,  *>  MODEL,  GRASL,  NPTS 
READ (5,*)  HI,  SMAX,  ANGLE 
DO  100  1=1, NPTS 

READ(  5,  *)  ZGR  <I),P(I)»T<I),RH(I),  WO  ( I ) , VIS ( I ) >  IHA(I) 

C  CONVERT  TEMPERATURE  TO  ABSOLUTE 

T< I )=T( I )+273.  15 

C  CHECK  AND  SET  DEFAULT  DATA 

IF  <P(I>  .  LT.  -500.0)  CALL  PRMOD(ZGR< I )+GRASL,  MODEL,  P < I > ) 

IF  (T(I)  .  LT.  -500.0)  CALL  TEMOD(ZGR( I )+GRASL, MODEL, T( I ) ) 

IF  (RH( I )  .  LT.  -500.0)  CALL  WHMOD ( ZGR ( I ) +GR ASL ,  MODEL ,  WH ( I ) ) 

IF  <WO(I)  .  LT.  -500.0)  CALL  WOMOD ( ZGR ( I ) +GRASL, MODEL,  WO ( I ) ) 

C  DETERMINE  WATER  DENSITY 

TA=273.  1 5/T  < I ) 

IF  <  RH  ( I )  .  LT.  -500.0)  GO  TO  50 
RHSAT=FCT  <TA) 

RH0=0.  01*RH< I) 

DN=1 .  0-(l.  O— RHO ) *RHSAT *R V*T < I ) /P < I ) 

WH ( I )=RHSAT*RHO/DN 
GO  TO  100 

C  CALCULATE  RELATIVE  HUMIDITY 

50  RHOSTR=P ( I ) / <  RV*T ( I ) ) 

RH< I ) =100.  0*< WH< I ) /FCT < TA) )*( ( RHOSTR-FCT ( TA ) ) / ( RHOSTR-WH ( I ) ) ) 

100  CONTINUE 
RETURN 


C**  PRMQD  **************■*****■*■****■*■»■•#•*■»••»•*•»•***■******•***■#■****•«■■«■*•«••»*■«•** 
C 

SUBROUTINE  PRMOD(Z,M,  P> 

C  PRESSURE  MODEL  DATA  ROUTINE 

DIMENSION  ZD(6) ,  PD(6,  6) 

C  ALTITUDE 


DATA 

*  ZD<1)  ,  ZD (2) 

*  O.  0  ,  1.  0  , 

PRESSURE  MODEL 

DATA 

*  PD<  1,  1 )  ,  PD(  1, 2)  , 

*  1.013E+03,  9.  040E+02, 

PRESSURE  MODEL 

DATA 

*  PD<2,  1 )  ,  PD ( 2,  2)  , 

*  1.  013E+03,  9.  020E+02, 

PRESSURE  MODEL 

DATA 

*  PD<3» 1 )  .  PD(3, 2)  . 

*  1.018E+03,  8.  973E+02, 

PRESSURE  MODEL 

DATA 

*  PD ( 4,  1  )  i  PD<4,  2)  , 

*  1.010E+03,  8.  960E+02, 

PRESSURE  MODEL 

DATA 

*  PD<5,  1 )  ,  PD ( 5,  2)  , 

*  1.  013E+03,  8.  878E+02, 

PRESSURE  MODEL 

DATA 

*  PD < 6.  1 )  ,  PD<6,  2)  , 

*  1.013E+03,  8.  986E+02, 


ZD(3) 
2.  O 


PD (1,3)  , 

8.  050E+02, 


PD(2,3>  , 

8.  020E+02, 


PD(3i  3) 

7.  897E+02, 


PD ( 4,  3)  , 

7.  929E+02, 


PD(  5»  3)  , 

7.  775E+02, 


PD  (6,3)  , 

7.  950E+02, 


SEMI -LOG  INTERPOLATION 
DO  lO  1=2,6 
L=  I 

IF  (Z  .  LT.  ZD(L) )  GO  TO  12 
CONTINUE 

FAC=(Z-ZD(L-1 ) )/(ZD<L)-ZD(L-l ) ) 
P=PD(M, L-l )*<PD(M,  L ) /PD (M,  L-l ) ) **FAC 
RETURN 

RETURN 


ZD  <  4 ) 
3.  0 


PD  (1,4)  , 

7.  150E+02, 


PD(2, 4)  , 

7.  100E+02, 


PD  <  3.  4 )  , 

6.  938E+02, 


PD ( 4,  4) 
7.  OOOE+' 


PD  (5,4)  , 

6.  798E+02, 


PD  (6,  4)  , 

7.  012E+02. 


ZD  ( 5 ) 
4.  0 


PD  (1,5)  , 

6.  330E+02, 


PD(2, 5)  , 

6.  280E+02, 


PD ( 3,  5) 

6.  081E+02, 


PD  (4,5)  , 

6.  160E+02, 


PD  (  5,  5 )  , 

5.  932E+02, 


PD ( 6,  5)  , 

6.  166E+02, 


ZD  ( 6 ) 
5.  0 


PD(  1,  6 
5.  590E 


PD(2, 6 
5.  540E 


PD  ( 3,  6 
5.  313E 


PD  ( 4,  6 
5.  410E 


PD ( 5,  6) 
5.  1 58E 


PD ( 6,  6) 
5.  405E 


279 

C 

280 

c**- 

TEMOD  **************************************************************| 

281 

c 

282 

SUBROUTINE  TEMOD(Z,M, T> 

283 

c 

TEMPERATURE  MODEL  DATA  ROUTINE 

284 

DIMENSION  ZD<6) ,  TE(6»  6) 

- 

285 

c 

ALTITUDE 

286 

DATA 

287 

*  ZD ( 1 )  *  ZD < 2 )  ,  ZD(3) 

ZD<  4) 

ZD  ( 5  > 

ZD ( 6 )  / 

288 

#  0.  0  ,  1.0  ,  2.  0  t 

3.  0 

4.  0 

5.  0  / 

289 

c 

TEMPERATURE  MODEL  1 

290 

DATA 

291 

*  TE<1,  1)  -  TE<  1 , 2)  ,  TE ( 1 i 3)  , 

TE  < 1 . 4 )  . 

TE (1,5)  , 

TE ( 1 , 6 )  / 

292 

*  3.  000E+02.  2.  940E+02,  2.  880E+02, 

2.  840E+02, 

2.  770E+02, 

2.  700E+02/ 

293 

c 

TEMPERATURE  MODEL  2 

294 

DATA  . 

295 

*TE(2.  1)  ,  TE<2, 2)  ,  TE(2,3>  , 

TE<2. 4)  , 

TE  ( 2,  5 )  , 

TE  <  2,  6 )  / 

296 

*  2.  940E+02,  2.  900E+02,  2.  850E+02, 

2.  790E+02, 

2.  730E+02, 

2.  670E+02/ 

297 

c 

TEMPERATURE  MODEL  3 

298 

DATA 

299 

*  TE ( 3,  1 )  ,  TE<3, 2)  >  TE<3,  3)  , 

TE  (3.4)  . 

TE(3.  5)  , 

TE  (3.6)  / 

300 

*  2.  722E+02,  2.  687E+02,  2.  652E+02, 

2.  617E+02. 

2.  557E+02, 

2.  497E+02/ 

301 

c 

TEMPERATURE  MODEL  4 

302 

DATA 

303 

*  TE<4>  1 )  -  TE<4, 2)  ,  TE<4,  3)  , 

TE  (4,4)  , 

TE  (4,5)  , 

TE  (4,6)  / 

304 

*  2.  870E+02,  2.  820E+02,  2.  760E+02, 

2.  710E+02, 

2.  660E+02, 

2. 600E+02/ 

305 

c 

TEMPERATURE  MODEL  5 

306 

DATA 

307 

*TE<5,1)  ,  TE( 5*  2 )  ,  TE<5,3>  , 

TE  (  5,  4 )  , 

TE  (  5,  5 )  , 

TE  (5,6)  / 

308 

*  2.  S70E+02,  2.  590E+02,  2.  559E+02. 

2.  527E+02, 

2.  477E+02, 

2. 409E+02/ 

309 

c 

TEMPERATURE  MODEL  6 

310 

DATA 

311 

*  TE (6.1)  ,  TE<6»  2)  .  TE(6, 3)  , 

TE(6, 4)  , 

TE<6. 5)  , 

TE  (6,6)  / 

312 

*  2.  881E+02,  2.  816E+02.  2.  751E+02, 

2.  687E+02, 

2.  622E+02, 

2.  557E+02/ 

313 

c 

SEMI-LOG  INTERPOLATION 

314 

DO  10  1=2.  6 

315 

L=I 

316 

IF  (Z  .  LT.  ZD(L) )  GO  TO  12 

317 

10 

CONTINUE 

318 

12 

FAC=  <  Z-ZD ( L— 1 ) )/<ZD(L)— ZD( L— 1 ) ) 

319 

T=TE(M.  L-1)*<TE<M,  L)/TE(M,  L-l)  )**FAC 

320 

c 

RETURN 

321 

RETURN 

322 

END 

323 

C 

324 

c** 

WHMOD  ***********************************************************##h*| 

325 

C 

326 

SUBROUTINE  WHMOD ( Z, M,  WH ) 

327 

C 

WATER  DENSITY  MODEL  DATA  ROUTINE 

328 

DIMENSION  ZD (6 ) i WHD (6,6) 

329 

C 

ALTITUDE 

330 

DATA 

331 

*  ZD  (1)  ,  ZD  (2)  ,  ZD  (3) 

ZD  (4) 

ZD  (5) 

ZD  (6)  / 

332 

*  0.  0  ,  1.  0  ,  2.  0  » 

3.  0 

4.  0 

5.  0  / 

333 

C 

WATER  DENSITY  MODEL  1 

334 

DATA 

335 

*  WHD  (1,1)  ,  WHD  (1,2)  ,  WHD  (1,3)  , 

WHD  (1,4)  , 

WHD (1,5)  , 

WHD  (1,6)  / 

336 

*  1.900E+01,  1.  300E+01.  9.  300E+00, 

4.  700E+00, 

2.  200E+00, 

1.  500E+00/ 

337 

C 

WATER  DENSITY  MODEL  2 

338 

DATA 

339 

«  WHD (2,1)  ,  WHD (2, 2)  >  WHD (2, 3)  , 

WHD  (2,  4)  , 

WHD(2,  5)  , 

WHD (2,  6)  / 

340 

*  1.  400E+01,  9.  300E+00,  5.  900E+00, 

3.  300E+00, 

1.  900E+00, 

1. OOOE+OO/ 

341 

C 

WATER  DENSITY  MODEL  3 

342 

DATA 

343 

*  WHD  (3,1)  ,  WHD  (3,  2)  ,  WHD  (3,  3)  , 

WHD (3,  4)  , 

WHD  (3,  5)  , 

WHD (3, 6)  / 

344 

*  3.  500E+00,  2.  500E+00,  1.  800E+00, 

1 . 200E+00, 

6.  600E-01, 

3.  800E-01/ 

345 

c 

WATER  DENSITY  MODEL  4 

346 

DATA 

347 

*  WHD (4,1)  »  WHD (4,  2)  ,  WHD (4,  3)  , 

WHD (4, 4)  , 

WHD (4,  5)  , 

WHD (4, 6)  / 

348 

*  9.  100E+00,  6.  OOOE+OO,  4.  200E+00, 

2.  700E+00, 

1.  700E+00, 

1. OOOE+OO/ 

349 

c 

WATER  DENSITY  MODEL  5 

350 

DATA 

351 

*WHD(5,1)  ,  WHD(  5,  2)  ,  WHD(5,3>  , 

WHD(5,  4)  , 

WHD (5,  5)  , 

WHD  (5,6)  / 

352 

*  1.  200E+00,  1.  200E+00,  9.  400E-01, 

6.  800E-01, 

4.  100E-01, 

2.  OOOE-Ol/ 

353 

c 

WATER  DENSITY  MODEL  6 

354 

DATA 

355 

*  WHD (6,1)  ,  WHD (6, 2)  .  WHD (6, 3)  , 

WHD (6,  4)  , 

WHD (6,  5)  , 

WHD (6,  6)  / 

356 

*  5.  900E+00,  4.  200E+00,  2.  900E+00, 

1.  BOOE+OO, 

1.  100E+00, 

6.  400E-01/ 

357 

c 

SEMI -LOG  INTERPOLATION 

358 

DO  10  1=2,6 

359 

L=I 

360 

IF  (Z  .  LT.  ZD(L) )  GO  TO  12 

361 

10 

CONTINUE 

362 

12 

FAC=( Z— ZD(L— 1 ) )/(ZD(L)— ZD(L-l) ) 

363 

WH=WHD(M, L-l )*( WHD(M< L)/WHD(M,  L-l ) ) **FAC 

364 

C 

RETURN 

365 

RETURN 

366 

END 

87 


r 


n  n 


**  WOMOD  #***************************************-k-****w***-#-***#******- 

SUBROUTINE  WOMOD  (  Z,  M*  WO ) 

C  OZONE  DENSITY  MODEL  DATA  ROUTINE 

DIMENSION  ZD ( 6 ) ,  WOD (6,6) 

C  ALTITUDE 

DATA 

*  ZD  (1)  ,  ZD  (2)  ,  ZD  (3)  ,  ZD  (4)  ,  ZD  (5)  ,  ZD  (6) 

*  0.  0  .  1.  0  ,  2.  0  ,  3.  0  ,  4.  0  .  5.  0 

C  OZONE  DENSITY  MODEL  1 

DATA 

*  WOD  (1,1)  ,  WOD  (1,2)  ,  WOD  (1,3)  ,  WOD  (1,4)  ,  WOD  (1,5)  ,  WOD  (1,6) 

*  5.  600E-05,  5.  600E-05,  5.  400E-05,  5.  100E-05,  4.  700E-05,  4.  500E-0 

C  OZONE  DENSITY  MODEL  2 

DATA 

*  WOD (2,  1 )  ,  WOD (2,  2)  ,  W0D(2,  3)  ,  W0D(2,4)  ,  W0D(2,5)  ,  W0D(2,6) 

*  6.  000E-05,  6.  000E-05,  6.  000E-05,  6.  200E-05,  6.  400E-05,  6.  600E-0! 

C  OZONE  DENSITY  MODEL  3 

DATA 

*  WOD (3,1)  ,  WOD (3, 2)  ,  WOD (3,  3)  ,  WOD (3. 4)  ,  WOD (3,  5)  ,  WOD (3, 6) 

*  6.  OOOE-05,  5.  400E-05,  4.  900E-05,  4.  900E-05,  4.  900E-05,  5.  SOOE-O! 

C  OZONE  DENSITY  MODEL  4 

DATA 

*  WOD (4,1)  ,  WOD (4, 2)  ,  WOD (4, 3)  ,  WOD (4, 4)  ,  WOD (4, 5)  ,  WOD (4, 6) 

*  4.  900E-05.  5.  400E-05,  5.  600E-05,  5.  B00E-05,  6.  000E-05,  6.  400E-0 

C  OZONE  DENSITY  MODEL  5 

DATA 

*  WOD (5,1)  ,  WOD  (5,  2)  ,  WOD (5,  3)  ,  WOD  (5,  4)  ,  WOD  (5,  5)  ,  WOD  (5,  6) 

*  4.  100E-05,  4.  100E-05,  4.  100E-05,  4.  300E-05,  4.  500E-05,  4  700E-0 

C  OZONE  DENSITY  MODEL  6 

DATA 

*  WOD (6,1)  ,  WOD (6, 2)  ,  WOD (6, 3)  ,  WOD (6, 4)  ,  WOD (6, 5)  ,  WOD (6, 6) 

*  5.  400E-05,  5.  400E-05,  5.  400E-05,  5.  000E-05,  4.  600E-05,  4  600E-0 

C  SEMI -LOG  INTERPOLATION 

DO  10  I»2. 6 
L=  I 

IF  (Z  .  LT.  ZD(L) )  GO  TO  12 
10  CONTINUE 

12  FAC=(Z-ZD(L-1 > ) / ( ZD ( L ) -ZD( L-l ) ) 

WO=WOD(M, L-l )*(WOD(M,  L)/WOD(M,  L-l ) )**FAC 
C  RETURN 

RETURN 


411 

C 

412 

c** 

413 

414 

C 

415 

416 

417 

418 

419 

C 

420 

421 

c 

422 

423 

424 

425 

c 

426 

100 

427 

c 

428 

110 

429 

430 

c 

431 

432 

c 

433 

434 

c 

435 

436 

c 

437 

438 

c 

439 

440 

441 

c 

442 

443 

444 

445 

446 

447 

448 

449 

c 

450 

c 

451 

452 

150 

453 

454 

455 

4  56 

c 

457 

458 

c 

459 

460 

c 

461 

462 

c 

463 

464 

c 

EQUABS  *****************  ****■«■•#  •*•&*■**•#■  *■***■#•■***■***■■»••»•*••»■***■*■*#*■■*■»•#■■«■•«■*•** 

SUBROUTINE  EQUABS<S,  IHAZE,  W, RHZ, TZ> 

CALCULATES  COMPONENT  AMOUNTS 

COMMON  /BLK1 /  ZGR ( 20 ) ,  P < 20 ) , T < 20 ) , RH < 20 ) , WH < 20 ) , WO ( 20 ) , VIS ( 20 ) 
COMMON  /BLK2/  HI  ,  ANGLE  ,  SMAX 

COMMON  /BLK3/  IHA<20),  NPTS 
DIMENSION  W<10>,  VSB<9) 

DETERMINE  ALTITUDE 
Z=H1+S*C0S(ANGLE*3.  14159/180.0) 

DETERMINE  DATA  ARRAY  LOCATION 
DO  100  1=2,  NPTS 
L=I 

IF  (Z  .  LT.  ZGR  <L ) )  GO  TO  110 
CONTINUE 

DETERMINE  FACTOR 

FAC=(Z-ZGR(L-1 ) ) / ( ZGR ( L ) -ZGR ( L-l ) ) 

PRESSURE  AT  Z 

PZ=P(L-1 )*<P<L)/P(L-1 ) )**FAC 
TEMPERATURE  AT  Z 
TZ=T ( L— 1)*(T<L)/T( L— 1 ) )**FAC 
WATER  DENSITY  AT  Z 
WHZ=WH(L-1 )*<WH<L)/WH<L-1 ) )**FAC 
OZONE  DENSITY  AT  Z 
WOZ=WO  <  L— 1 ) * ( WO ( L ) /WO ( L— 1 ) )**FAC 
RELATIVE  HUMIDITY 
RHZ=RH<L-1)*(RH(L)/RH(L-1 ) >**FAC 
SET  AEROSOL 
IHAZE=IHA<L-1 ) 

VISZ=VIS<L-1) 

EQUIVALENT  AEROSOL  AMOUNT 
W  ( 7 )  =0.  0000 

IF  ( IHAZE  .  EQ.  0)  GO  TO  150 
DATA 

*  VSB <  1 ) ,  VSB <  2  > ,  VSB <  3  > ,  VSB ( 4 ) ,  VSB ( 5 ) ,  VSB ( 6  > , VSB (7), VSB ( 8 ) , VSB ( 9 ) / 

*  23.  O  ,  5.  0  ,  23.  0  ,  5.  0  ,  5.  0  ,  50.  O  ,  23.  O  ,  0.  2  ,0.5  / 

IF  (VISZ  .  LT.  -500.0)  VISZ=VSB ( IHAZE ) 

W ( 7 ) =3.  912023/VISZ 

PA=PZ/ 1013.  0 
TA=273.  15/TZ 

SET  EQUIVALENT  AMOUNTS 
D=0  1*WHZ 
X=PA*TA 

PT=PA*SQRT<TA) 

EQUIVALENT  H20  AMOUNT 
W< 1 )=D*PT**0.  9 

EQUIVALENT  C02  N20  AND  ETC.  AMOUNT 
W(2)=X*PT**0.  75 

EQUIVALENT  N2  AMOUNT 
W ( 4 )  =0.  8*PT*X 

PPW=4.  560E-05*D*TZ 


19 


TS 1=296. O/TZ 

EQUIVALENT  H20  CONTINUUM  AMOUNT  (10  UM) 
W(5)=D*PPW*EXP<6.  08*<TS1-1.  0) >+0. 002*D* ( PA-PPW > 
EQUIVALENT  H20  CONTINUUM  AMOUNT  (4  UM) 

W  <  9 ) =D* ( PPW+O.  12* ( PA— PPW ) ) *EXP  <4.  56*(TS1-1.  0) ) 
EQUIVALENT  MOLECULAR  SCATTERING  AMOUNT 

W(6)=X 

EQUIVALENT  OZONE  AMOUNT  (U.  V.  REGION) 

W  <  8 ) =46.  6667*W0Z 

EQUIVALENT  OZONE  AMOUNT  (I.R.  REGION) 
W(3)=W(8)*PT**0.  4 
RETURN 


478  C 

479  C**  FCT  #•■*■#■***■*■** *#*** *****«•**■**#■  **■«■***  *#■«•**  w***##* ***** ■#■*<•** *#** *■*■*•** 

480  C 

481  FUNCTION  FCT ( A ) 

482  C  FUNCTION 

483  FCT=A*EXP < 18.  9766-14.  9595*A-2.  43882*A*A) 

484  RETURN 

485  END 


91 


486  C 

487  C**  PLAN  ********************w*********************-********#**-a'#****** 

488  C 

489  FUNCTION  PLAN(T,  IV) 

490  C  PLANCK  BLACKBODY  FUNCTION 

491  XLAM=1  0E+04/IV 

492  C1=0.  59544E+08 

493  C2=14388.  000 

494  AUG=C2/<XLAM*T> 

495  PLAN=2.  0*C 1 / ( < XLAM**5 ) * ( EXP ( AUG ) - 1 .  0) ) 

496  RETURN 

497  END 


I? 


C**  C 1DTA  *****#****#**#**********************************«***********'»' 
C 

SUBROUTINE  ClDTAtCl,  IV) 

C  DETERMINE  WATER  VAPOR  SPECTRAL  COEFFICIENT 

COMMON  /SPEC1/  C1ARR (2580) , C2ARR ( 1375 ),C3ARR( 540).  C4ARR< 133) 
COMMON  /SPEC2/  C7ARR<280,  5),  C8ARRU02.  2)  »C9ARR(131) 

1*0 

Cl-O.  000 

IF  (IV  .  GE.  350  .AND.  IV  .  LE.  9195)  I-(  IV-350) /5+1 

IF  (IV  .  GE.  9875  .AND.  IV  .  LE.  12795)  I«( IV-9875)/5+1771 

IF  (IV  .  GE.  13400  .AND.  IV  .  LE.  14520)  I*(  I V- 13400)/ 5+2356 

IF  (I  .EG.  0)  RETURN 
Cl»10.  000**C 1 ARR ( I ) 

RETURN 


529  C 

530  C*«  C3DTA  *#********#*#**###*#*#******#******##*#*#*####*******##****## 

531  C 

532  SUBROUTINE  C3DTA(C3. IV) 

533  C  OZONE  I.  R.  SPECTRAL  COEFFICIENT 

534  COMMON  /SPEC1/  C 1ARR< 2580 ), C2ARR< 1575).  C3ARR< 540),  C4ARR< 133) 

535  COMMON  /SPEC2/  C7ARR< 280, 5),  C8ARR(102,2)  ,C9ARR<131) 

536  I»0 

537  C3*0.  000 

538  IF  <IV  .  GE.  575  .AND.  IV  .  LE.  3270)  I-(  IV-575)/5+l 

539  IF  <1  .  EQ.  0)  RETURN 

540  C3-10.  000**C3ARR  < I ) 

541  RETURN 

542  END 


C**  C40TA  *****#*#*#*****#*#***#**♦**###**####*******##*##*#*#******** 
C 

SUBROUTINE  C4&TA<C4,  IV > 

C  NITROGEN  CONTINUUM  SPECTRAL  COEFFICIENT 

COMMON  /SPEC1/  C1ARR(2980)< C2ARR( 1979)*  C3ARR(940)«  C4ARR( 133) 
COMMON  /SPECS/  C7ARR<380» 9),  C8ARR<102»2>  * C9ARR (131) 

I«0 

C4-0.  OOO 

IF  (IV  .  GE.  2080  .AN©.  IV  .  LE.  2740)  I-(  IV-2080) /5+1 

IF  (I  .  CQ.  O)  RETURN 

C4»C4ARR(I) 

RETURN 


"557  C 

558  C**  C9DTA  #****♦#***#*###*#*****#************##**####************** 

559  C 

560  SUBROUTINE  C5DTA(C5,  IV) 

561  C  WATER  VAPER  CONTINUUM  (10UM)  SPECTRAL  COEFFICIENT 

562  C5-0.  0000 

563  IF  (IV  .  GE.  670  .AND.  IV  .  LE.  1350)  CO  TO  50 

564  RETURN 

565  50  V-IV 

566  C5-4.  18+5578.  0*EXP<-7.  87E-03*V) 

567  RETURN 

568  END 
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%  *'*  . V.** 


v 


O  <0 


C**  C60TA  *********#♦♦*###**###*###**#***#*#**####****#*****#******* 
C 

SUBROUTINE  C6DTA<C6,  IV) 

C  MOLECULAR  SCATTERING  SPECTRAL  COEFFICIENT 

C6-0.  000 

IF  <IV  .  GE.  4000  )  00  TO  50 
RETURN 
50  V-IV 

CA«V**4/<9.  26799E+18—1.  07123E+09*V**2) 

RETURN 


981 

982 

983 

984 
989 

986 

987 

988 

989 

990 

991 

992 

993 

994 
999 

996 

997 

998 

999 
600 
601 
602 

603 

604 
609 
606 

607 

608 

609 

610 
611 
612 

613 

614 
619 
616 

617 

618 

619 

620 
621 
622 

623 

624 
629 
626 

627 

628 

629 

630 

631 

632 

633 

634 


C 

C**  C7DTA  ************************************************************* 
C 

SUBROUTINE  C7DTA(C7E,  C7A, IV. RH. .IHAZE) 

C  DETERMINE  AEROSOL  SPECTRAL  COEFFICIENT 

COMMON  /SPEC1/  C1ARR (2980) »  C2ARR ( 1979) . C3ARR ( 940) .  C4ARR  < 133) 
COMMON  /SPEC2/  C7ARR (280.  9) .  C8ARRU02.  2)  .C9ARRU31) 

DIMENSION  RHZ0NE(4) 

DATA 

«  RHZONE ( 1 ) .  RHZ0NE(2) .  RHZ0NE(3).  RHZ0NE(4)/ 

*  0.  000.  70.  000,  80.  000.  99.  000/ 

C  CHECK  FOR  NO  AEROSOL 

C7A-0.  000 
C7E«0.  000 

IF  (IHAZE  .  EQ.  0)  RETURN 
C  DETERMINE  DEPENDENT  VARIABLES 

V*IV 

Y-l.  0E+04/V 
X-100.  0“RH 

C  CHECK  FOR  FOGS 

IF  (IHAZE  .  EQ.  8  .OR.  IHAZE  .EG.  9)  GO  TO  lOOO 
C  DETERMINE  DATAFILE  COLUMN 

DO  100  K-2.  4 
J-K 

IF  (RH  .LT.  RHZONE ( J )  )  GO  TO  lio 
100  CONTINUE 
110  Xl-100.  O-RHZONE(J-l) 

X2-100.  O-RHZONE(J) 

C  ALLOW  FOR  WAVELENGTH  IN  FIRST  COLUMN 

J-J+l 

C  DETERMINE  EXTINCTION  COEFFICIENT 

C  DETERMINE  TYPE  OF  AEROSOL 

IF  (IHAZE  .  EQ.  1  .OR.  IHAZE  .  EQ.  2)  ISTRT-2 

IF  (IHAZE  .EQ.  3  .OR.  IHAZE  .  EQ.  4)  ISTRT=162 

IF  (IHAZE  .  EQ.  9)  ISTRT-82 
IEND-ISTRT+38 

C  SET  UP  WAVELENGTH  INTERPOLATION 

DO  200  K-ISTRT,  I END 
I-K 

IF  (Y  .  LT.  C7ARR ( 1 ,  1 ) )  GO  TO  210 
200  CONTINUE 

C  LOO  ON  RH  AND  LINEAR  ON  WAVELENGTH 

210  Y1-C7ARR(I— 1, 1) 

Y2-C7ARR(I, 1) 

ZA1*C7ARR  < I— 1» J-l ) 

ZA2*C7ARR(I-1,  J) 

ZB1»C7ARR(I,  J-l ) 

ZB2-C7ARR(I,  J) 

FAC A-ALOQ ( Z A2/ Z A 1 ) / ALOG ( X2/ X 1 ) 

FACB*ALOO ( ZB2/ZB1 ) /ALOG ( X2/X1 ) 

ZA-ZA1* ( X/X 1 ) **FACA 
ZB-ZB1*(X/X1 )**FACB 
FAC* ( Y— Y 1 ) / ( Y2-Y 1 ) 

C7E-FAC# ( ZB-ZA ) *ZA 


A  • 


DETERMINE  ABSORPTION  COEFFICIENT 

1*1+40 

LOG  ON  RH  AND  LINEAR  ON  WAVELENGTH 
Y1*C7ARR  <1-1* 1) 

Y2*C7ARR ( I* 1) 

ZA1*C7ARR ( I-li  J-l) 

ZA2*C7ARR< I-l»  J) 

ZB1*C7ARR ( I» J-l) 

ZB2*C7ARR(I.  J) 

FACA*ALOG( ZA2/ZA1 > /Al_(30< X2/X 1 ) 

FACB-ALOO ( ZB2/ZB 1 ) /ALOG ( X2/X 1 ) 
ZA*ZA1*'X/X1  >**FACA 
ZB-ZB1*(X/X1>**FACB 
FAC*  <  Y-Y 1 ) / <  Y2— Y 1 ) 

C7A-FAC*  < ZB-ZA  >  +ZA 
RETURN 

FOGS 

DETERMINE  EXTINCTION  COEFFICIENT 
1000  IF  ( I  HAZE  .EQ.  8)  J*2 
IF  ( IHAZE  .EQ.  9)  J-4 
ISTRT-242 
IEND-ISTRT+38 
DO  1100  K-ISTRT.  IEND 
I*K 

IF  (Y  .  LT.  C7ARR ( I >  1 ) )  GO  TO  1110 
1100  CONTINUE 
lllO  Y1*C7ARR< I— 1» 1) 

Y2“C7ARR  ( I «  1 ) 

Z1*C7ARR( I-li  J) 

Z2*C7ARR(  I.  J) 

FAC— < Y— Y1 ) / ( Y2— Y1 ) 

C7E— FAC* ( Z2— Z 1 ) +2 1 

DETERMINE  ABSORPTION  COEFFICIENT 
IF  (IHAZE  .EQ.  8)  J-3 
IF  (IHAZE  .EQ.  9)  J-3 
Yl— C7ARR ( I— 1» 1) 

Y2—C7ARR ( I « 1) 

Z1*C7ARR( 1-1.  J) 

Z2— C7ARR( I.  J) 

FAC*(Y— Yl )/(Y2— Yl ) 

C7A-FAC* ( Z2-Z 1 ) +Z 1 
RETURN 


C#*  C8DTA  **#*#******************■)»•■■»■-. 

C 

SUBROUTINE  C8DTA(C8,  IV) 

C  OZONE  U.  V.  SPECTRAL  CG£~ 

COMMON  /SPEC1/  C 1 ARR  <  2580 ) »  C2.-V 
COMMON  /SPEC2/  C7ARR(280,  5),  C3.- 
C8-0.  000 

IF  (IV  .  GE.  13000  .AND.  IV  .  LE. 

IF  (IV  .  GE.  27500  .AND.  IV  .  LE. 
RETURN 

C  13000  CM-1  TO  24000  CI'I- 1 

100  V*IV 

DO  120  1*2 .  56 
N-I 

IF  (V  .  LT.  C8ARR  (N»  1 )  )  GO  TO  12' 
120  CONTINUE 
N*56 

130  FAC*  ( C8ARR  ( N— 1 »  2>-C8ARR<N,  2)  )  /  (  C 
C8»FAC*(V-C8ARR(N-1,  1 ) 1+C8ARR  L  - 
RETURN 

C  27500  CM-1  TO  50000  CM-1 

200  V*IV 

DO  220  1*58, 102 
N*I 

IF  (V  .  LT.  C8ARR  (N,  1 )  )  GO  TO  22  j 
220  CONTINUE 
N-102 

230  FAC- ( C8ARR ( N-l ,  2 ) -C8ARR ( N*  2 )  )  /  (  C 7 
CB-FAC* ( V-C8ARR ( N-l ,  1 ) > +C3ARR < f  i- 
RETURN 


•-.-***■*«*****#**********##*###*# 


:  zn  t 

f/5). C3ARR ( 540 ) > C4ARR ( 133) 
.102,2)  i C9ARR ( 131 ) 

■.  >0 )  GO  TO  100 

::oo>  go  to  200 


:  R ( N- 1 ,  1 ) -C8ARR ( N,  1  )  ) 

>: ) 


•7?  (N-l,  1 )-CBARR(N,  1 )  ) 


c 

SUBROUTINE  C9tfTA<C9, IV) 

C  NATER  VAPOR  CONTINUUM  <4UM>  SPECTRAL  COEFFICIENT 

COMMON  /SPEC.l/  C 1ARR  (2980)  /  C5ARR  ( 1J573)  *  C3ARR  ( 540 ) .  C4ARR  ( 133 ) 
COMMON  /SPEC2/  C7ARR 1280,  3 ) ,  C8ARR(102. 2)  ,  C 9ARR(131) 

1*0 

C9*0.  000 

IF  (IV  .  GE.  2390  .AND.  IV  .  LE.  3000)  I*< IV-2330) /5+1 
IF  (I  .£Q.  0)  RETURN 
C9*C9ARR< I ) 

RETURN 

END 


•  ’  V  •  '  N  ' 


